Tree structures for the current fluctuations in the exclusion process 
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We consider the asymmetric simple exclusion process on a ring, with an arbitrary asymmetry 
between the hopping rates of the particles. Using a functional formulation of the Bethe equations of 
. the model, we derive exact expressions for all the cumulants of the current in the stationary state. 

' These expressions involve tree structures with composite nodes. In the thermodynamic limit, three 

\ regimes can be observed for the current fluctuations depending on how the asymmetry scales with 

. the size of the system. 
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I. INTRODUCTION 



■ The one-dimensional asymmetric simple exclusion process (ASEP) is a stochastic model featuring classical hard-core 
particles hopping between neighboring sites of a one-dimensional lattice, with an asymmetry between the hopping 
& • rates forcing a current of particles to flow through the system. The ASEP thus belongs to the large class of driven 
diffusive systems which have been playing an important role in the understanding of of out of equilibrium 

1 , systems. In particular, it is a special case of the Katz-Lebowitz-Spohn model [J, [f| describing a lattice gas with 
particles subject to nearest-neighbor interactions and driven by an external field. Since it is one of the simplest 
examples of an interacting particles model with a non equilibrium steady state, the ASEP has been studied much in 
the past (Ml. 

Various boundary conditions have been considered in the study of the ASEP. Open boundary conditions [Tl| - tl3 | 
have been used to model the coupling of the system to reservoirs of particles. If the particles are interpreted as quanta 
of energy, the ASEP models heat transport between two heat baths at different temperatures [l5| . The ASEP has 
also been considered on an infinite line and on a finite ring fl9r{23l ]. Even though the model defined on the 

infinite line can be seen as the infinite system size limit of the periodic model, these two models behave differently 
\Q , in their respective stationary state. Indeed, the stationary state corresponds to taking the infinite time limit, which 
£C) ■ does not always commute with the infinite system size limit. In this paper, we will consider only the periodic model 
on a ring. 

The ASEP can be related to several other models of statistical physics, in particular models of a fluctuating 
' interface growing through the deposition and evaporation of particles. The deposition of a particle in the growth 
. model corresponds in the ASEP to the move of a particle in the forward direction, while the evaporation of a particle 
^ ' corresponds in the ASEP to the move of a particle in the backward direction. Thus, the fluctuations of the integrated 
• i-h . current in the ASEP are related in the corresponding growth model to the fluctuations of the height of the growing 
' interface. The symmetric exclusion process (SSEP), for which the forward and backward hopping rates are equal, 
H I corresponds to a fluctuating interface with the same deposition and evaporation rates. It is described at large scales 
by the Edwards- Wilkinson equation (24|. On the other hand, the asymmetric exclusion process, with different forward 
and backward hopping rates, corresponds to a growing interface described at large scales by the Kardar-Parisi-Zhang 



equation [25| . The SSEP can then be seen as discrete version of a system evolving by the Edwards- Wilkinson equation, 
while the ASEP is a discrete version of the Kardar-Parisi-Zhang equation. 

The fluctuations of the current of particles is an important quantity for the exclusion process. In the case of the 
periodic system on a ring, these fluctuations have been calculated using mainly two methods: the matrix Ansatz [2fj| . 
first introduced for the ASEP in for the calculation of the stationary measure of the open system, and the Bethe 
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Ansatz, which relies on the underlying quantum-integrability of the model. The diffusion constant, which characterizes 
the average deviation of the current from its mean value, has initially been calculated using the matrix Ansatz: first 
in [22| for the the totally asymmetric model (TASEP), for which the particles hop only in the forward direction, and 
then in [28| for the more general case of the partially asymmetric model. Higher cumulants of the current have been 
obtained subsequently using the Bethe Ansatz. In [29(, the thermodynamic limit of all the cumulants was obtained 
for the symmetric model. Before that, in (30l. l3l| . finite size expressions have been calculated for all the cumulants 
of the current in the totally asymmetric case. In |32| . this result was generalized in the thermodynamic limit to the 
partially asymmetric case with non-vanishing asymmetry. However, this result did not allow to study precisely the 
transition between the symmetric and the totally asymmetric regimes, as this transition occurs in a scaling where the 
asymmetry vanishes in the thermodynamic limit. It was thus needed to find finite size expressions for the cumulants 
of the current as this would allow to study all the possible scalings for the asymmetry. A first step in this direction 
was to recover the exact expression for the partially asymmetric diffusion constant from the Bethe Ansatz, using a 
functional formulation of the Bethe equations [33j . This calculation was then extended to the third cumulant of the 
current in [34||, and all the cumulants of the current were obtained in [35| for the weakly asymmetric model (WASEP), 
for which the asymmetry scales as the inverse of the size of the system. 

In this article, we obtain finite size expressions for all the cumulants of the current in the model with arbitrary 
asymmetry, generalizing the known results for the three first cumulants, as well as the result for all the cumulants in 
the totally asymmetric model. We write the cumulants of the current in terms of sums over tree sets. These expressions 
for the cumulants, that we conjectured earlier in (36{, are proved here using functional Bethe Ansatz equations. From 
these exact expressions, we study the thermodynamic limit of the cumulants of the current in various scalings for the 
asymmetry. We observe in particular that the Edwards- Wilkinson and the Kardar-Parisi-Zhang regimes arc separated 
by an intermediate regime corresponding to an asymmetry large with respect to the inverse of the size of the system 
but small with respect to the square root of the inverse of the size of the system. 

The article is organized as follows. In section [Hi we summarize the main results of this paper for the cumulants 
of the current. In section we recall the Bethe Ansatz for the exclusion process and prove that the solution 
of the functional Bethe equation corresponding to the stationary state is solution of a simpler functional equation. 
Then, in section IIV1 we show that this simpler functional equation gives rise to tree structures. In section [V] we 
express the solution of the functional Bethe equation in terms of these tree structures. In section I VII we write a 
parametric expression for the generating function of the cumulants of the current. Finally, in section fVIH we reduce 
this parametric expression to an explicit expression for the cumulants of the current involving forest structures. A 
few technical calculations are relegated to the appendices. 



We consider the asymmetric simple exclusion process on a ring of size L with n particles hopping locally both one 
site forward (with rate p, i.e. the probability of hopping is pdt in any infinitesimal time interval of length dt) and 
backward (with rate q). By the exclusion rule, the particles are only allowed to hop if the destination site is empty. 
This dynamics is represented in fig. Q]for a configuration with three particles on sixteen sites. In all this article, we 
will use the notation 



for the ratio of the hopping rates. The symmetric model then corresponds to x = 1, while the totally asymmetric 
model corresponds to x — 0. In the following, we will call 1 — x the asymmetry between the hopping rates. 

In this section, we will recall known results about the current fluctuations in the stationary state of the asymmetric 
exclusion process on a ring. We will then summarize the results that will be obtained in the rest of the article. 



We call Y t the total distance covered by all the particles between time and time t. Each time a particle moves 
forward, Y t increases by 1, while each time a particle moves backward, Y t decreases by 1. The quantity Y t is thus 
the integrated current of all the particles. In the long time limit, the probability to observe a value j of the current 
different from its mean value vanishes exponentially fast (l5| as 



II. CURRENT FLUCTUATIONS IN THE ASYMMETRIC EXCLUSION PROCESS 



A. Cumulants of the stationary state current 




(2) 



FIG. 1: Transition rates for the asymmetric exclusion process on a ring. The particles move one site forward with rate p and 
one site backward with rate q if the destination site is empty. 



The function G is called the large deviations function of the current. Introducing a fugacity 7 associated to Y tl we 
define the Legendre transform E(j) of G(j) by 

£7(7) = max[j7 - G(j)} ■ (3) 
The long time behavior of the mean value of e lYt is then given in terms of £7 (7) by the relation 



(4) 



The previous equation means that £7 (7) is the exponential generating function of the cumulants of the current in the 
stationary state. If we write its perturbative expansion near 7 = as 

„, , , . D(x) , E*(x) q Ea(x) a , , 

£(7) = J(z)7 + -jrl 2 + ^p7 3 + -^p7 4 + ■ ■ • , (5) 

then the derivatives of £(7) at 7 = give 

J(x) = lim - — (6) 

i-voo t 

D(x) . lim ghg . lim MJ*> (7) 

t— >-oo £ *— >oo t 

e 3{ x) = lim m^mm+wi = lim <^-w) 8 ) (8) 

E , x) = Um - 4(y t )(r t 3 ) - 3(y t 2 ) 2 + i2(y t ) 2 (r t 2 ) - e(y t ) 4 

i->oo t 

= lim ((y t -(y t )) 4 )-3((y t ~(y- t )) 2 ) 2 

t— S-oo f 

Here, J(x) is the mean value of the current, D(x) the diffusion constant, and E%{x) and E±(x) respectively the third 
and the fourth cumulant of the current. 

The mean value of the current J(x) can be calculated using the stationary measure of the exclusion process on a 
ring || P s tat(C) = 1/ (^) for any configuration C of the system. This leads to 

J M = {l - x) < L - n ) , (io) 

p L — 1 

The higher cumulants of the current are more difficult to obtain, as they involve correlation functions at different times 
and not just equal-time correlation functions. In particular, they can not be obtained knowing only the stationary 
measure. The diffusion constant D(x) was first calculated by Derrida and Mallick in [28[ using an extension of the 
matrix Ansatz for the stationary state. It was then recovered by Bethe Ansatz in [33| . It is given by 

D(x) _ 2(1 - x)L ^ 1 + x k ( n ^ k ) ( n ^ fc ) 



p L-l ^ { lj 
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The third cumulant of the current was then obtained in [34[ using again Bethe Ansatz. The exact result for E$(x) 
can be written as 



(L — l)E 3 (x) 1 ^ V"^ (-2 i • • , -2\ (n+i) (n+j) in-i-i) 



L \l L W L 
n I 

L \( L \( L 



(12) 



L \ ( L \ ( L \f L 



, 2 s L L +i )L-i)L L +j )L-o) i + 



1 - a;l«l 1 - scIj'I 



where all the (1 + x'°')/(l — a;' ') must be replaced by an arbitrary constant A which cancels from the result. 



B. Combinatorial formula for the cumulants of the current 



In this article, we will obtain finite size expressions f| 125[) for all the cumulants of the current, generalizing the 
previous expressions pT|) and (|12p of the diffusion constant and the third cumulant. After solving the functional 
Bethe equation at all order in the fugacity 7, we will prove that the fc-th cumulant Ek(x) has the following structure: 



Ek (x) 1 — x ( L 
2 



V 



L - 1 



fc-i 



W(h) 
Sf(h) 



(13) 



In this expression, Hk-i is a particular set of forests (a forest being a set of trees) that will be defined in section 
IVII Al The first sets Hk are given by 
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In equation ([T3|) . the rational number Sf(h) is a symmetry factor associated to the forest ft,, while is equal to 

k — 1 nested sums over the integers of product of binomial coefficients and factors of the form 



1-xH 



(18) 
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For example, for the fourth cumulant of the current, we will find 



(L - l)E 4 (x) 
p(l - x)L 3 




' 1 +2:1*1 1 + x\ j \ 1 + xl fe l 1 + 
l-a;H 1-xW 1 - 



x (i 2 +ij +f +jk + k 2 ) 






) 



(19) 



In the totally asymmetric limit x = 0, this expression leads to (using some binomial coefficient formulas, see e.g. the 
appendix A of [34|) 



In the symmetric exclusion process (x = 1), the forward and backward hopping rates are equal and the system 
on a ring reaches an equilibrium stationary state in the long time limit. On the contrary, the partially asymmetric 
exclusion process (i ^ 1) reaches in the long time limit a non equilibrium stationary state. These two system thus 
belong to different universality classes. It is known [lj[ that at large scales, the dynamics of the symmetric model can 
be described by the Edwards- Wilkinson (EW) equation [24j . while the dynamics of the partially asymmetric model 
can be described by the Kardar-Parisi-Zhang (KPZ) equation [25j . 

The crossover between the EW and the KPZ regimes can be studied more precisely by looking at models for which 
the asymmetry 1 — x scales in the large system size limit as 1 — x ~ \/L r for a positive real number r. It can be 
expected that for large values of r, the model lies in the EW regime, while for small values of r, the model lies in 
the KPZ regime. Finding the values of r for which the system belongs to the EW or KPZ regimes is thus a natural 
question. 

A heuristic argument (35j indicates that both values r = 1 and r = 1/2 correspond to a natural separation between 
a weakly and a strongly asymmetric model. This tends to show that there are in fact three distinct regimes for the 
current fluctuations: a regime 1 — x -C 1/L corresponding to the EW equation, a regime 1 — x ^> 1/ \[L corresponding 
to the KPZ equation, and an intermediate regime (I) for which the current fluctuations are neither described by the 
EW equation nor by the KPZ equation. The scaling 1 — x ~ 1/L is usually called the weakly asymmetric scaling. In 
the following, we will call the scaling 1 — x ~ 1/vL the strongly asymmetric scaling. 

The weakly asymmetric scaling has received much attention recently. In this scaling, a phase transition has been 
observed by Derrida and Bodineau [37L |38| using a hydrodynamical approach to the current fluctuations based on the 
"macroscopic fluctuations theory" of Bertini, De Sole, Gabrielli, Jona-Lasinio and Landim [39[. There is a critical 




(20) 



which agrees with the exact solution [30( of Derrida and Lebowitz. 



C. Three regimes for the current fluctuations 
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value v c of the asymmetry such that for v < v c , the generating function of the cumulants of the current E(j) is 
quadratic when L — > oo, while for v > i/ c , the function £'(7) is not quadratic even when L — > 00. This phase 
transition corresponds to a change in the density profile adopted by the system when a current different from the 
mean value of the current J(x) is forced to flow through the system. The phase transition can also be seen from a 
non analyticity in the generating function of the cumulants of the current in the weakly asymmetric scaling [35 1. 

The existence of the three regimes EW, I and KPZ can be justified by taking the thermodynamic limit L, n — > 00 
with particle density p = n/L fixed in the the exact formulas for the cumulants of the current. For the diffusion 
constant, it was shown in [28[ that the exact formula (jTTJ) becomes for large systems 



f ~^(1- x)p 3/2 (1 - p) 3/2 L 3/2 if ^ « 1 



f~43p(l 



P)LJ ° 



2 -u* 

U "tanh($ M ) 



2p(l - P)L 



if 1 



2$ 



if 1 — x < 



1 



KPZ regime 
contains TASEP 



EW and I regimes 
contains SSEP 



(21) 



We observe on these expressions that the behavior of the system is separated in two distinct regimes, depending on 
whether 1 — x is large or small with respect to 1/ y/~L. The regime 1/VL< 1 — x is the KPZ regime, which contains 
in particular the totally asymmetric model. The regime 1 — x <C 1/ yL correspond to the reunion of the EW and I 
regimes, and contains the symmetric model. We note that the diffusion constant keeps the same value in the EW and 
I regimes. It becomes larger in the KPZ regime, where it depends on the asymmetry. 

For the third cumulant of the current, it was shown in [34| that in the thermodynamic limit 



Ms. 
p 



s E3 



-(1 



8tt 
3^3 



3tt 
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3 2 (1 



f ~2$/i 3 ($)p 3 / 2 (l-p) 3 / 2 L 5 / 2 



2 (1 



Here, h^i^) is given by the double integral 



p) 2 L 3 if -4= < 1 - x (KPZ regime) 



p)*(l-Z B p(l-p))L 



if 1 



2$ 



if-^<l-a;<C^ (I regime) 



(22) 



if 1 



x ~ T 



if 1 - x <T 4j (EW regime) 



ha(*) = - 



2tt 
375 



-'o 



o?m dv 



(u + v 2 )e 



1\ n -u*-v* 



(u + uv + v )e 



2\ n —u —uv- 



tanh($u) tanh(<f>w) 



(23) 



Contrary to the case of the diffusion constant, the three regimes EW, I and KPZ give different expressions for the 
third cumulant of the current. In particular, the third cumulant is positive in the EW regime and becomes negative 
in the I and KPZ regimes. The presence of three regimes is also observed on the higher cumulants of the current. In 
[1^, [35j], all the cumulants were obtained in the EW regime, as well as in the weakly asymmetric scaling 1 — x ~ 1/L. 
Previously [30l - l32l |. all the cumulants were obtained in the KPZ regime. From the exact expression (|125j) . we will 
obtain in section IVIIDI the expression of all the cumulants in the strongly asymmetric scaling 1 — x ~ 1/ \f~L. The 
only regime left is the intermediate regime. 

We observe that in this regime, the third cumulant of the current is simply given by the limit v — > 00 of the 
expression in the weakly asymmetric scaling, which agrees with the limit $ — > of the expression in the strongly 
asymmetric scaling [34[ . We checked that the limits v — > 00 and $ — > give also the same expression in the case of 
the fourth cumulant. We conjecture that this is true for all the cumulants of the current, which allows to obtain the 
expression in the intermediate regime by taking the limit v — > 00 of the expression in the weakly asymmetric scaling 
(as it is much more difficult to take the limit $ — >• of the expression in the strongly asymmetric scaling). 

We can finally write the following expressions in the various scaling limits for the higher cumulants of the current 
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E k (k > 3): 



^ ~ (1 - z)/i fc (°°)P (fe+1)/2 (l - p )(fe+i)/2 L 3(fe-i)/2 

^ ~ 2$/ lfe ($) j0 fc / 2 (l - p )k/2 L (3k-4)/2 
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(fc even) 
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if 1 _ x „ 2* 

y/p(l-p)L 

if \ « 1 - a; « ^ (I) 
ifl-x-f 

if 1 - x < i (EW) 



(24) 



The Bj are the Bernoulli numbers, and the functions h k are given in equation (|134j) by A: — 1 nested integrals. There 
are two distinct expressions for E k in the EW regime depending on the parity of k. This is related to the fact that 
for the symmetric model, the symmetry between the forward and the backward directions on the ring forces the odd 
cumulants of the current to be equal to zero, unlike the even cumulants. 



III. BETHE ANSATZ FOR THE ASYMMETRIC EXCLUSION PROCESS 

In this section, we write the functional Bethc equation whose polynomial solution Q gives access to the cumulants of 
the stationary state current in the exclusion process. Then, we recall the construction due to Pronko and Stroganov of 
a polynomial P solution of a "dual" functional Bethe equation. Combining the polynomials P and Q into a function 
w, we finally obtain a closed equation for the unknown quantity w. This last equation will be solved perturbatively 
to all order in the fugacity 7 in the next sections. 



A. Reminder of the functional formulation of the Bethe equations 

The generating function of the cumulants of the current E(^) is the eigenvalue with largest real part of a deformation 
M(7) of the Markov matrix of the system [3(J [Hj]. The matrix M{^) is related to the Hamiltonian of the XXZ spin 
chain defined in terms of the usual spin 1/2 operators by 



H 



XXZ 



\ E (Sfs^ x + S^S^ + Afl$*>flfM ) , (25) 



with the "twisted" boundary condition (25 (±) = S {x) ± iS^) 



c( + ) _ ( ^ Y q(+) _ ( ^ X \ o(-) CrOO _ cr(*) 



and with a parameter A given by 



A = V " 3 > 1 ■ (27) 

More precisely, Hxxz is related to M("f) through the similarity transformation (see e.g. [l(| for the case 7 = 0; the 
generalization to nonzero 7 is straightforward) 

Hxxz ^M( 7 )-L^l. (28) 

Since the Hamiltonian Hxxz is known to be intcgrable, the matrix M(j) is diagonalizablc using the Bethe Ansatz. We 
note that E(j), which is the largest eigenvalue of M(^), also gives the ground state of the (non-hermitian) quantum 
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Hamiltonian i?xxz- 

Using the Bethe Ansatz, each eigenvalue E of M (7) can be expressed in terms of a solution of the Bethe equations 
of the system 



L 



l^jM = _TT«zfw (29) 

1 - xy t J AJ. ajy. _ y. 



as 



The Bethe equations (|29[) form a set of n coupled polynomial equations for the quantities y\ , . . . , j/„ . Different 
solutions of this set of equations give through ([50]) different eigenvalues of M(7). Introducing the polynomial 

n 

Q(t) = H(t- Vj ), (31) 
i=i 

the Bethe equations can be rewritten [33| as the functional equation 

Q(t)R(t) = e i7 (l - t) L Q{xt) + x n (l - xt) L Q(t/x) , (32) 

where the two unknown polynomials Q and R are of respective degrees n and L. This functional equation is the scalar 
version of Baxter's equation [io| . As the Bethe equations (j2T?)) , it has several solutions corresponding to different 
cigenstatcs of the matrix M(7). The solution corresponding to the largest eigenvalue of M(pj) is the only one such 
that 

Q(t)=t" + 0( 7 ) . (33) 

This solution also verifies [331 



z"Q(l/z) 



0(1) 

The corresponding eigenvalue E(j) is given by 



e™ 7 . (34) 



E{1) - -(!-*)# log f^g^) (35) 



p v ' dt °\ Q(t) 



|t=i 



The functional equation p2p has been used in [3J] to calculate the three first cumulants of the current. In the present 
article, we will show that another functional equation, equivalent to (|32p . is more suitable to calculate systematically 
all the cumulants of the current. 



B. Functional equation "beyond the equator" 



Given a polynomial Q of degree n < L/2 solution of the functional Bethe equation (|32|) . it is possible to construct 
a polynomial P of degree L — n such that 



P(0) 



(36) 



and 



(1 - .T™e- L7 )Q(0)(1 - t) L = Q{t)P(t/x) - x n e-^Q(t/x)P(t) 



(37) 



The construction of the polynomial P, due to Pronko and Stroganov |4lj |. is explained in appendix|A] As the functional 
Bethe equation (|3^|) . equation (|3T|) has several solutions corresponding to several eigenvalues of the deformed Markov 
matrix M (7). In the following, we will calculate the solution Q(t) corresponding to the largest eigenvalue £'(7) by 
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solving (|37j) instead of (1321) . From equations (f33|) . (l36l) and (|37|) . this solution is characterized by Q(t) = t n + O (7) 
and 

P(t) = l + 0( 7 ). (38) 

We note that using (|37|) to replace Q(xt) and Q(t/x) in equation (|3"21) by expressions in terms of both polynomials P 
and Q gives an equation in which all the Q(t) cancel: 

P(t)R(t) = x n (l - t) L P(xt) + e L7 (l - xt) L P{t/x) . (39) 

From this last equation, we observe that (P, x L ~ n e~ L ~ t R) is in fact solution of the initial functional Bethe equation 
(|32l) with n replaced by L — n and 7 replaced by log x — 7. 



C. Resolution of the functional Bethe equation 

The functional equation (|37|) depends on two unknown polynomials P and Q, which makes it difficult to solve 
directly. In this section, we will gather the polynomials P and Q corresponding to the stationary state into a unique 
function w which preserves all the information from both P and Q. We will then rewrite the functional equation (|37j) 
as a closed equation (|59[) for the quantity w. This equation, which was not known before, is the key to a systematic 
calculation of all the cumulants of the stationary state current in the exclusion process on a ring. 

The characterization ([33]) and (|38]l of the polynomials P and Q corresponding to the stationary state makes it 
natural to solve the functional equation (|37| perturbatively near 7 = 0. We write the perturbative expansion of Q(t) 
and P(t) near 7 = as 

Q(t) =t n + Q 1 (t.) 1 + Q 2 (t) 1 2 + ... (40) 
P(t) = 1 + P a (t) 7 + P 2 (t) 7 2 + • ■ ■ (41) 

Since the term of higher degree of Q{t) is equal to t n , all the Qk(t) are polynomials in t of degree n — 1. Similarly, 
since P(0) = 1, all the Pfc(0) are equal to zero. In the perturbative expansion near 7 = 0, the quantity 



\ t n J t n \ t n 2t 2n 

has thus only strictly negative powers in t, while the quantity 



(P 2 (t)-^)7 2 + 



has only strictly positive powers in t. The fact that both \og(Q(t)/t n ) and log(P(t)) can be expressed as a sum of 
powers of t in such a way that no power of t appears in both quantities, which is a consequence of the perturbative 
expansion near 7 = 0, will be crucial in the following. We write 

Q{t)\ _ ^ xH j 



j<0 



x ] V 



log(P(t))=^[/3] JT — -, (45) 



j>0 



where the [a]j and the [jS\j are formal series in 7. The sums over j cover respectively the strictly negative integers 
and the strictly positive integers. At each order in 7, only a finite number of j give a nonzero contribution to the 
sums. From the coefficients [a]j and [(3]j we also define 

a(t) = (46) 

j<0 

/3(t)^J2^t 3 ■ (47) 
i>o 

In the previous expression, a(t) has only strictly negative powers in t: at each order in 7, it is a polynomial in 1/t 
without constant term. Similarly, (3(t) has only strictly positive powers in t: at each order in 7, it is a polynomial 
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in t without constant term. We now introduce the operator X acting on any series with both positive and negative 
powers in t 



*(*) = XH-** 



(48) 



as 



with the convention 



X[u(t)] = J2[u} j t j - 



1 + 



- x\i\ 



1 + arl°l 



A 



(49) 



(50) 



l-a?l°l 

for an arbitrary constant A. Because of the absolute value in its definition, the operator X acts differently on quantities 
with only negative powers in t and on quantities with only positive powers in t. Using (|44[) and the definition (|46[) of 
ait), we find in particular the following relations: 

a(t) X[a(t)] 



log 



x n Q(t/x) 



2 2 
a(t) X[a{t)\ 



t n J 2 

Similarly, using (|4"5"j) , and the definition (|4T|) of Pit), we obtain 



log(P(*)) = -f + «1 



log (P(i/a;)) 



Pit) , X[/3(f)] 



(51) 
(52) 

(53) 
(54) 



2 2 

The four previous equations allow us to express the rhs of the functional equation (|37[) only in terms of ait), Pit) and 
the operator X: all the dependency on the asymmetry x has been absorbed in these three quantities. The functional 
equation ([57)1 now becomes 

(1 — t) L c(t)-g(t) x[ Q (t)-3(t)] 



(1 



r.n e -Lj\ 



Q(o)- 



-L 7 + 



(55) 



We observe that the previous equation depends on a and /3 only through the difference a{t) — Pit). Thus, we define 

ait) L 7 ^t) 



wit) 



Using the notations 



and 



hit) = 



t n 



2 IP 2 — TP 2 ) 

C = — { - — — ^ Q(0) 



and recalling that X[l] = A, equation (15511 finally becomes 



w(t) = arcsinh (Chit) 



,X[w(t)] 



(56) 



(57) 



(58) 



(59) 



This closed functional equation depends only on one unknown quantity wit), unlike the functional equation (|37[) which 
involved two unknown polynomials P and Q. This is the key to the systematic calculation of all the cumulants of the 
current that will be performed in the following. 

We note that in the initial functional equation (|3"T)) , the fact that P and Q are polynomials was crucial to constrain 
the set of the solutions: for a given polynomial P of degree L — n, there is in general no polynomial Q of degree n 
such that P and Q form a solution of (|3"T|) . In equation (|5"9"|) however, we do not need to use the fact that a and f3 
are constructed using polynomials anymore: we already used it to prove that a(t) has only negative powers in t and 
Pit) has only positive powers in t. 
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D. Remaining steps for the calculation of the generating function E(y) 



Expanding equation (|59[) as a formal series in the parameter C, we note that it admits a unique power series 
w(t) = O (C) as a solution. At first orders in C, this solution can be obtained directly from the expansion in C of 
(j5"5f . Up to order 5 in C, we have 



){t) =hC + hX[h]C 2 
T 



-hX\h} 2 + hX[hX[h}] - -h 3 
2 6 



1 



hX[h] 3 + hX[h]X[hX[h}] + -hX[hX[h] 2 ] + hX[hX[hX[h]]] - -h 3 X[h] - -hX[h 6 ] G 



— hXlh} 4 
24 L J 



1 



1 



-hX[hX[h]] 2 + -hX[h} 2 X[hX[h}} + -hX[h]X[hX[hY 



2 1 L " 2 
hX[h}X[hX[hX[h}}} + hX[hX[h]X[hX[h}]] 



6 



-hX[hX[h} 3 } 



-hX[hX[hX[h] 2 ]] + hX[hX[hX[hX[h]}]] 



-hi 3 X[h] 2 - Ui 3 X[hX[h]] - ^hX[h]X[h 3 } - ^hX[h 3 X[h}} - ^hX[hX[h 3 }} + ^h 5 



(60) 



where we wrote h for h(t) to lighten the notations. In section ITVl we will perform a systematic calculation of w(t) to 
all order in the parameter C. This calculation will involve tree structures. Indeed, we note on the previous equation 
that each term is of the form a rational number times a function depending on t through h(t) and some nesting of 
X operators. It is convenient to write such a term as a tree representing the nesting of the X operators. The edges 
of the tree are labeled by X and the nodes are labeled by h to the power an odd integer. For example, we have the 
correspondence 



\hX[h]X[h 3 ] 



\hX[h] 2 X[hX[h]] o 




(61) 



Then, once w(t) is known, a(t) can be obtained by extracting its negative powers in t using (|56[) . This is done in 
section [V] But, from (|4"4"]) and we have 

, , , fx n Q(t/x)\ 

Using and ([55)1 . we will thus obtain a parametric expression for E(^) in section fVTl equation will give £^(7) 
in terms of the parameter C, while equation (|34[) will give 7 in terms of C . The parameter C is of order 7 from (|58[) 
and (|33[) . It means that in the parametric expression for £'(7), the parameter C can be eliminated order by order in 
7. This elimination is done systematically in section fVIII using the Lagrange inversion formula. It will finally give an 
explicit expression for all the cumulants of the current in terms of forest structures. 



IV. SOLUTION OF THE CLOSED FUNCTIONAL EQUATION IN TERMS OF TREE STRUCTURES 

In this section, we will show that the closed functional equation for the quantity w(t) defined in (f5B|) can be 
used to obtain an explicit expression of w(t) involving tree structures. We will first define these tree structures and 
then show how they appear from the functional equation 1)59 p . 
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A. Rooted trees with composite nodes 

We call "elementary node" an object that will be represented as a dot •. The elementary node is the brick 
with which we will build the tree structures necessary to express the fluctuations of the current in the asymmetric 
exclusion process. We call gencrically " composite node" an object containing an odd number of elementary nodes. In 
this section, we will consider only composite nodes without internal structure on the elementary nodes they contain. 
The size |c| of a composite node c is the number of elementary nodes it contains. 

We define the set C of composite nodes without internal structure. A composite node c G C will be represented by 
a set of dots corresponding to the elementary nodes it contains, surrounded by a closed line. We have 

c = |®' ®' @) •••} ' (63) 

where we have only drawn the four composite nodes of C up to size 7. 

We now build trees whose nodes are composite nodes. The size of a tree g, which will be noted \g\, is defined 
as the sum of the sizes of the composite nodes of g. It is also the total number of elementary nodes contained in 
the composite nodes of g. We define the set Q° of the rooted trees with nodes belonging to C. For a tree g £ G°, 
one elementary node is called the root (or elementary root) and will be represented by a small circle O instead of a 
black dot. The composite node containing the elementary root is called the composite root. For a tree g element of 
Q°, we call e(g) the set of the elementary nodes of g and c(g) the set of the composite nodes of g. For a composite 
node c € c(g), we define the number v c of the composite nodes of g which are neighbors of c (that is, the number of 
composite nodes linked to c by an edge). 

For each positive integer r, we also define the set Q° of the trees in Q° of size r. The five first sets Q° arc 
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The first sets Q° can be explicitly constructed using a computer. For r between 1 and 16, we find that the number of 
trees in Q° is given by 



r 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 
card£° 1 1 3 6 15 36 94 245 663 1815 5062 14269 40706 117103 339673 991834 ' 

In the following, we will see that the trees of Q° correspond to the good tree structures to express w(t). 

B. to (t) as a sum over rooted trees 

The expansion of the arcsinh in equation (|59[) gives 



(65) 



r=l ^ 
(r odd) 



(66) 



with rll — r x (r — 2) x (r — 4) x . . . x 3 x 1 for any positive odd integer r. We write the expansion of w(t) as a formal 
series in C 



(t) = X> fe (t)C fe , (67) 



it) _ 

fc=i 

and we expand e rX l w wi in powers of C . This expansion generates tree structures, as we have seen in (|60|) and (|61[) 
up to order 5 in C. In appendix [Bj we prove from equation (|66j) that 

This expression involves the set of trees defined previously. The factor x(,9) can be computed by induction on a 
tree g using 



dg) = Hty x\x(g')] 



g subtree of g 
connected to the 
composite root 



(69) 



if the tree g is made of a composite root of size r to which the trees g' are attached (the product is taken to be equal 
to 1 if the tree g is made of a single composite node). The symmetry factor S (g) is given by 

(number of equivalent^ 
elemSr/oo^of J 




\ elemental v rooi oi u / /i~rr\\ 

x ' (70) 



s o(9) \ ±* M 3 |c|! P(g) 



where P(g) is the number of permutations of the composite nodes of g leaving g invariant (except for the choice of 
the composite root). We recall that c(g) is the set of the composite nodes of g, and v c the number of composite nodes 
neighbors of c. 

We note that S a (g) depends on the position of the root of g only from the number of equivalent choices of the 
elementary root. This will allow us to replace the sum over the set of rooted trees in Wk(t) by a sum over a set of 
unrooted trees, whose number is smaller than the number of rooted trees. 



C. Unrooted trees with composite nodes 



We call Q the set of unrooted trees whose nodes are composite. The set Q can be obtained from the set Q° of rooted 
trees by replacing the elementary roots of the trees by ordinary elementary nodes. Several rooted trees from Q° will 
then correspond to the same unrooted tree from Q. 

For a tree g element of Q, we call again e(g) the set of the elementary nodes of g and c(<?) the set of the composite 
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nodes of g. For a composite node c £ c(g), we also define the number v c of the composite nodes of g which are 
neighbors of c. We call Q r the set of unrooted trees of size r, that is, the set of unrooted trees containing r elementary 
nodes. The five first sets Q r are 



a = {®} & = {®-®} & = {®-®-®.M 



04 = <®-®-®-®,<§ 



0b H®-®-®-®-® 




» 



(71) 



In the following, we will use a simplified notation for the trees of Q. All the composite nodes which contain only 
one elementary node will simply be represented as a dot •, while the other composite nodes will be represented by 
their size surrounded by a circle (e.g. (3), (5), . . . ). We emphasize that with these simplified notations, a dot • will 
represent a composite node of size 1 and not an elementary node. Using the new notations, the six first sets Q r are 




06 = 



(72) 



The first sets r can be explicitly constructed using a computer. For r between 1 and 16, we find that the number of 
trees in Q r is given by 



r 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 
card0,. 1 1 2 3 6 12 25 55 126 304 745 1893 4893 12916 34562 93844 



(73) 



D. w (t) as a sum over unrooted trees 
We will now rewrite the expression (|68[) of Wk (t) using a sum over g £ Gk instead of the sum over g £ G^- We have 



w k (t) 



geg k 



g has the same 
tree structure as g 



x(g') 
So(g') 



(74) 



where the trees g' are all the trees of 0£ that can be obtained from g by choosing an elementary root. Calling gj the 
tree of G% obtained from g by choosing for elementary root the j-th elementary node of g (for some arbitrary order 
on the elementary nodes of g), we can also write 



1 



X(9j 



(^) ^ ^ ^ ^ /numbe 

g eg k j = l I choices of the u °\y3J 
\ elementary root of gjj 



(75) 
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9 eg 

S(g) 
P(g) 

geg 

S(g) 
P(g) 



geg 

S(g) 
P(g) 



geg 

8(g) 
P(g) 



-18 



6 -6 



1 2 



•— •— «^ • o • (| ) • » •— (3)— • @ 



24 
24 



-4 200/3 
2 1 



120 



120 



^ • • • 

-6 
1 



dHD 



-2 
1 



-12 -4 40/3 72 

2 6 12 



TABLE I: Symmetry factors of the trees of g up to size 6, as given by equation (|76 



From the expression (|70|) of the rooted symmetry factor S {gj), we define an unrooted symmetry factor S(g) associated 
to the unrooted tree g by 



(76) 



We recall that c(g) is the set of the composite nodes of g, v c the number of composite nodes of g neighbors of c. and 
P(g) the number of permutations of the composite nodes of g that leave g invariant. With this unrooted symmetry 
factor, we finally obtain 




(77) 



The symmetry factors of the unrooted trees up to size 6 are given in table Q] 
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V. EXACT PARAMETRIC SOLUTION OF THE FUNCTIONAL BETHE EQUATION 



In this section, we write explicitly the expansion of w(t) in powers of C and t, and give a parametric expression of 
the polynomial Q(t). In order to do that, it is useful to add an internal tree structure to the composite nodes of the 
trees. 



A. Composite nodes with internal tree structure 

In section HVAl we defined the set C of composite nodes without internal structure. We now define the set C of the 
composite nodes with an internal unrooted tree structure linking all the elementary nodes contained in the composite 
node. A composite node c € C will be represented by a set of dots • representing the elementary nodes it contains, 
linked by edges representing the tree structure on the elementary nodes, and surrounded by a closed line. For the five 
composite nodes of C up to size 5, we have 




We recall that the size of a composite node (i.e. the number of elementary nodes it contains) is required to be an odd 
integer. 

We now define the set Q of the rooted trees with elementary nodes belonging to C, oriented edges labeled by i±, 
i2, • ■ ■ , the edges between composite nodes being identified to edges between elementary nodes belonging to different 
composite nodes. The labels of the edges of g £ Q are chosen all different. For g £ G, one node is called the root (or 
elementary root) of g and is represented by a small circle O instead of a dot •. Again, the composite node containing 
the elementary root will be called the composite root. We have 




where we represented all the trees of Q up to size 3. We recall that the size of a tree is equal to the number of 
elementary nodes it contains. As before, we introduce the set Q r made of the trees of Q of size r. Equivalently, the 
set Q r is the set of the trees of Q with r — 1 edges. The four first sets Q r are (again, we do not draw circles around 
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the composite nodes of size 1 to simplify the notations) 




The . . . represent the trees of Q r which can be obtained from the ones drawn by permuting the labels of the edges or 
reversing the directions of the edges. 

As in the case of the trees of G, we define for a tree g £ Q the set e(g) of the elementary nodes of g and the set c(g) 
of the composite nodes of g. For example, the four first trees of G4 drawn in ((80)) have four composite nodes (all of 
size 1), while the seven following trees have one composite node of size 1 and one composite node of size 3. 

In the following, we will make a distinction between two types of edges in the trees: "inner edges" , which link 
elementary nodes belonging to the same composite node, and "outer edges" which link elementary nodes belonging 
to distinct composite nodes. We will also call "inner labels" the labels of the inner edges and "outer labels" the labels 
of the outer edges. The set of the outer edges of a tree g will be noted o(g). For example, the four first trees of Q4 
drawn in (|80p have three outer edges and no inner edge, while the seven following trees have one outer edge and two 
inner edges. 



B. Function denned on the trees 

We will now define functions acting on the sets Q or Q, from which the cumulants of the current will be expressed 
in the following. 

Let g be a tree of Q. For an outer edge o € o{g), the subtree of g beginning at o is the tree made of the edge o 
and of all the edges which can be reached from o by moving on the edges of g (independently of the direction they 
are pointing to) away from the elementary root of g. We then call m{o) the linear combination (with coefficients +1 
and —1) of the indices ij which label the edges (inner and outer) of the subtree of g beginning at o, with coefficients 
chosen in the following way: the labels of the edges pointing in the direction of the root have a coefficient —1 while 
the labels of the edges pointing away from the root have a coefficient +1. We illustrate this definitions on the example 
of the tree 




This tree g has four outer edges 03, 04, 07 and o$, labeled respectively by 13, i-j and i%. They are such that 
771(03) = h, 771(04) = U - h + i& + ii — h, ?7i(o 7 ) = i? and 77i(o 8 ) = — is- 

We also want to associate a linear combination of the indices ij to each elementary node. For an elementary node 

e e e(g) of a tree g G Q (e may be the elementary root of g), we call £(e) the sum of the labels of the edges (inner 
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or outer, linking e to another elementary node) pointing to e, minus the sum of the labels of the edges pointing away 
from e. For the nine elementary nodes e € e(g) of the tree g drawn in (|8ip . the £(e) are 7 3 , i\ — 7 3 , — i\ — 72, 7 2 — *4j 
74 + k - k,_~k - k, «7, «6 + «8 and^-7 8 . 

For g £ Q, we define a tree g* G by attaching all the edges (inner and outer) of g to the elementary root. The 
edges keep their labels, their direction toward the elementary root (an edge pointing in the direction of the root in g 
will still point in the direction of the root in g* , and conversely), and their inner or outer nature. All the composite 
nodes of g* are thus of size 1 except possibly the composite root. For example, if we consider the tree g drawn in 
(f8"Tj) . we have 




(82) 



For this tree g* , if we call 03, 04, 07 and 0% the outer edges labeled by 23, 74, 77 and 7s, we have 777(03) = 73, 
771(04) = H, 777(07) = k and m(os) = —is- For the nine elementary nodes e G e(g*), we also have the following 
values for the £(e): ii, 12, 23, 74, —75, 76, 17, —is (e elementary node different from the elementary root of g*) and 
— ii — 72 — 73 — 24 + 75 — jg — 77 + is (e elementary root of g*). 

We will now define functions mapping a tree of Q or £7 to a number. For g € G, and 93 and 77 two arbitrary numerical 
functions, we define 



u VlV ig) = ( £ ^( £ ( e )) n 

eee(g) I \e£e(g) 



(83) 



For the trees g and g* represented respectively in (|8Tj) and (|82|) , we have for example 

U v . v {g) = [<p{h) + <p(h - h) + f{-ii - i%) + <p(i2 - u) + ¥>(u + k - k) 

+ <p(-i5 ~ ij) + fih) + f{k + k) + (f(-ia)] 
x v{h)v(h - k)v(-h - k)il(k - u)v(k + k~ k)v(~k - ^vfaWk + k)v(~k) 
U v , v {g*) = [<p(ii) + <p(i2) + ¥>(k) + v(U) + f(-k) + <p(h) + ¥>(k) 

+ ip(-k) + <p(-ii - 7 2 - 2 3 - 2 4 + 7 5 - 7 6 - 7 7 + i 8 )] 
x v(ii)v(i2)v(k)v(k)v(-k)v(k)v(k)v(-k)v(-ii ~h~ k - U + k - k - k + k) 



(84) 



(85) 



By definition, U v ^ v (g) depends only on the tree structure (with directed edges) of the elementary nodes of g, and 
is independent of the position of the elementary root and of the composite nodes in g. For an arbitrary numerical 
function £, we also define 



(86) 



For a tree g without any outer edge (i.e. a tree formed of only one composite node), we take the product in V^(g) to 
be equal to 1. For the trees g and g* represented in (|5T|) and (|82|) . we have 

V € {g) = £(2 3 )£(i 4 - 75 + 2 6 + k - k)i{k)i(-k) (87) 
V s (g*) = ZikWiKikKi-k) ■ (88) 

Unlike U VtTI (g) , V£ (g) depends on the position of the composite root, of the composite nodes of g, and on the orientation 
of the edges, but not on the internal tree structure of the elementary nodes inside the composite nodes (except for 
the direction of the edges). 

We now choose a map 9 from Q to Q preserving the tree structure of the composite nodes; the map 9 applied on 
g £ Q roots g, adds an internal tree structure on the elementary nodes of each composite node, transforms the edges 
between composite nodes into edges between elementary nodes, and adds a direction and a label to each edge. Thus, 
the map 9 changes a tree from Q into a tree from Q . 

For example, the first tree of drawn in (|72p can be mapped by 9 to the first or the second tree of Q4 drawn in 
([8TJ|) , and also to every tree that can be obtained from it by changing the orientation of the edges and any permutation 
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of the labels. In the same way, the second tree of G4 can be sent by 9 to the third or fourth tree of Q4. Finally, the 
third tree of Q4 can be mapped by 9 to the fifth up to eleventh tree of G4 drawn in (f80l) . 
We now define for g G Q r 



(89) 



with the convention W£'^(9) = 99(0)77(0) . This definition of W<2 (<?) is equivalent (see appendix [D|) to the definition 



w™ c G?) = £•••£ £W%nn(%)) 



(90) 



It can be proved (see appendix [Ej that the function W£'* does not depend on the choice of the map 9 if the function 
£ is even, which will be the case in all the situations we will consider in the next sections. We emphasize that if the 
function 77 has a finite support, the sums over the indices £1, . . . , £ r _i in the definition of W^(g) are finite sums. 
In the following, the generating function of the cumulants of the current will be expressed in terms of the function 



C. Parametric expression of the polynomial Q 



In appendix IFl we expand the expression ([77} of w(t) in powers of t. Then, keeping only the negative powers in 
t, we obtain an expression for the polynomial Q(t) in terms of the function WR* defined previously. We define the 
functions ipi, 77 and £a by 



T}(z) = 



(n4 


-*) ( 


(L- 


-0 


L » 




n-\-Z/ 
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u 









L — n—z\ 



if z ^ 



(91) 
(92) 
(93) 



Then, the polynomial Q solution of the functional Bethe equation (|3"2"j) corresponding to the stationary state (|3"3"j) is 
given by 



In 



x n Q{t/x) 

Q(t) 



00 00 



(94) 



where B, which depends on the arbitrary parameter A unlike Q{t), is equal to 



e-T* ) Q(0) 



(95) 



This expression for Q(t) is precisely the one that we conjectured in |36j . where it was checked up to order 5 in 7 
and 1 — t for all the systems up to size L = 12 by using the systematic pcrturbativc solution of the functional Bethe 
equation (|32|) presented in (34|. 

We note that divergences appear in equation (|94[) in the terms such that L < I, because of the denominator of <pi(z). 
An analytical continuation for complex L shows that these divergences cancel in fact in the limit where L becomes 
an integer as soon as <pi{z) is multiplied by rj(z), which is systematically the case in W^f x (g). The terms such that 
L < I in (|94j) must then be understood by taking the analytic continuation for complex L and by taking the limit L 
integer. 

We emphasize that the function £^ is even. This ensures that the previous expressions do not depend on the choice 
of the map 9 which enters the definition of W^f x {g). 
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VI. PARAMETRIC EXPRESSION OF THE CURRENT FLUCTUATIONS 

In this section, we use the solution of the functional Bethe equation derived previously to give a parametric 
expression of the generating function of the cumulants of the current E(y) . Then, we show that this expression allows 
to recover the known exact result in the totally asymmetric model |30ll. and the result in the partially asymmetric 
model with non- vanishing asymmetry in the thermodynamic limit [32J . 



A. Parametric form of the generating function Ei^y) 

Taking t = 1 in equation (|94|) gives, using ([34]) . an expression for 7 in terms of B: 



1 v^V v w ^ x{9) 



(96) 



The function tpo is the affine function 



<Po(z) 



L 



(97) 



In the definition of U ViV (g) for a tree g <E G, a summation over all the elementary nodes e of g of the cp(£(e)) is 
performed. But, by definition of £(e), the summation over e of the £(e) is equal to zero. Indeed, for each edge of 
<7, there is only one node to which the edge is pointing toward and one node from which it is pointing away, which 
implies that each index ij appears in one of the £(e) with a coefficient +1 and in another one with a coefficient — 1. 
The linear term of ip then does not contribute to U v 



nor to Wg't* 



Wc obtain 




(98) 



We have written z h> 1 for the function equal to 1 everywhere. Taking now the derivative in t = 1 of equation (|94| . 
the expression (|35[) for the generating function of the cumulants of the current implies that 



Bin) 



= -2{l-x)Y J 



B 



E 



2 J S(g) 



(99) 



The function ipi is a polynomial of degree 2: 



n(L — n) (L — 2n)z z 2 



L(L-l) L(L-l) L(L-l) 



(100) 



As previously for ip , the linear term in epi does not contribute to W^f x (g), whereas the constant term gives again 
We finally find 




(101) 



Here, z 1— > z 2 is the square function, and the mean value of the current J is given by (flU)) . Because of the definition 
of W applied to the only tree of Qi, the term k = 1 does not contribute to the sum in the previous equation. 

Equations (|98|) and (jlOip give a parametric expression of the generating function of the cumulants of the current: 
the parameter B can be eliminated between both equations to give an expansion of ^(7) in powers of 7 which does 
not involve B. We will see in section IVIII that the elimination of the parameter B can in fact be performed in a 
systematic way to all orders in 7. 
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B. Totally asymmetric limit 

We note that the parametric expression (|98[) . (|101|) of the generating function of the cumulants of the current E(j) 
looks very much like the one obtained by Derrida and Lebowitz [30T ] for the totally asymmetric model. We will now 
prove that the result of [3(| can be recovered from and (|101j) . In the limit x 0, the function £>, becomes 



(£0^0 ( z ) 



A if* = (1Q2) 
1 ifz^O V 



We choose the arbitrary parameter A to be equal to 1 so that £\(z) = 1 for each value of z. From (|86[) . it implies that 
V^ A (0(g)) is equal to 1 for any tree g € Qk- We also note that if the map 8 defined in section IV Bl is chosen so that all 
the edges in 8(g) point away from the elementary root, then 8(g)* is the same tree for all g € Qk (we recall that the 
map 8 transforms a tree of Q into a tree of Q). From ([83]) . and with the function 77 defined in (|92|) , we find 

U v ,r,(8(g)*) = [<p(h)+ ■ ■ ■ + f(ik-i) + <p(—ii - ... - ifc-i)) 

( t )•••( / ) x ( L ) 

x Vra+Zfc.i/ Vw-ti-...-» fc _ 1 ; (103) 



Thus, neither W^^* (g) nor W^TlZ (g) depend on the tree g anymore. They can be calculated by considering the 



zi—tz 

generating function (1 — ti) L ... (1 — tk) L ■ We obtain 



rkL\ 

W^ 1 (g) = k^L (104) 



= k (k-lML-n) G& 



fr},z\—*-l 



The parametric expression f[95|) , fjlQip of E(j) for the totally asymmetric model then rewrites 

^ = -Z^^T7I^ x 2.^) ( 106 ) 

£( 7 )-J 7 = n(L-n) ^ k(k-l)B k (g) / v 1 \ 

L(X-l)^(JfeL-l)2*-i • ( j 

In appendix [Cl we calculated a generating function for the trees of Qk with the symmetry factor S. In z = 1, equations 
(|Clj) and (|C2[) give the sum over the trees in the previous equation: 

Expressing the term J7 of E("f) using the expression of 7 in terms of B, we finally obtain 

00 k fkL\ 

fc=i (j 

m_ n(L-n) ™ g* g) 

^ fe=l U 

By absorbing into B the binomial coefficient to the power k, we recover the result of [30j . 



C. Thermodynamic limit with finite asymmetry 



For a system with particle density p = n/L fixed, it is possible to take the thermodynamic limit of the equations 
(I98[) and (|101[) giving £'(7) under a parametric form if the asymmetry 1 — x stays finite when L — > 00. This allows 
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to recover the result obtained by Lee and Kim [32J in this particular limit. The steps of the derivation are similar to 
the case of the totally asymmetric limit studied in the previous section. From the large L expansion 

(1 — n \ Z z 2 + (l-2p)3 
e " 2 P( i-p)i. , (111) 

and from the expression ([90]) of (fiO> the parametric expression (|98|) and (jlOip for £"(7) becomes for large L 

„ OO / s fc ^ z 2 W 

t—zE^J £ ^ (112) 

fc=l V 7 3 6S fc VJJ 

k W z ^ x{z) z2 (g) 

^ ^ fc=2 v 7 gSSfc Vay 

We used the fact that the terms in e z of rj(z) cancel in W^(g), which is similar to the cancellation of the linear term 
of if that we used in section IVl Al 

Because of the —z 2 /L in the exponentials, the indices ij in W contributing the most to the sums are of order y/L. 
Thus, we set ij = Uj^J Lp(l — p). The Riemann sums over the indices ij in W can now be rewritten as integrals. We 
note from equation (|93[) that if A is equal to 1 and 1 — x remains strictly between and 1 in the limit L — ¥ 00, we can 
replace £\{z) by 1. It implies that we can replace V^ x (9(g)) by 1 in the Riemann integrals, using the definition ([90)1 of 
W . We also note that if the map 9 is chosen so that all the edges of 9(g) point away from the elementary root, then 
Utp,r)(0(g)*) does not depend on the tree g anymore. It means that we can use the generating function of the trees 
(|C2j) as in the case of the totally asymmetric limit. We obtain 



fc— 1 _ 



ifc-1 



t_i ™ J -co J -co 

- — («! + ... +«J_ 1 +«l«2+«l«3 + ---+«fe-2«fc-l) 



(114) 



fe=2 

2(tt? + . . . + U 2_ i + WiU2 + UiU3 + _ _ + Ufc _ 2Ufc _ 1 ) e -Cu;+.»+«J_ 1 +t. I u 8 +UiUs+...+u*- a u k _ 1 ) _ ( U5) 



We can now perform the integration over the indices uj. Setting finally D = —B^/2irp(l — p)L, we find the result of 
Lee and Kim 1321 at first order in L: 



7' 



1 ^ (-£>) 



2ttp(1 - p)L 3 ^ /fc 3 / 2 



p - g V 2ttL 3 \ fc 3 / 2 fc 5 / 2 



In this limit, the dependency in the asymmetry in £"(7) reduces to the global factor p — q = p(l — x). 



VII. EXPLICIT EXPRESSION FOR THE CUMULANTS OF THE CURRENT 



We wrote in the previous section a parametric expression for the generating function of the cumulants of the current 
E(j) of the partially asymmetric model. This expression involved sums over tree sets. In this section, we will give 
another expression for £"(7), equivalent to the previous one, as an explicit formal series in 7. It will involve sums over 
forest sets. This new expression for £(7) will give explicit expressions for the cumulants of the current. In particular, 
we will recover the known results for the three first cumulants. 



A. Forests with composite nodes 

We call "/oresi" an acyclic graph whose connected components are trees. For a forest h, we will note h the number 
of trees in h and \h\ the size of h, defined to be the sum of the sizes of all the trees it contains. The size of a forest is 
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thus equal to the sum of the sizes of all its composite nodes. It is also equal to the total number of elementary nodes 
it contains. A forest will be represented as a list of trees surrounded by square brackets [ ]. We will consider in the 
following two sets of forests. 

We first define the set Ti of the forests made of trees elements of Q of size larger or equal to 2. We call Tt r the 
subset of T-L made of the forests h such that \h\ — h = r. The four first sets Ti r are 



n 2 = {[• • •] 

n 3 = 

U A = 







(118) 



, [CD-m] , [•-CD-*] , [(s 



»,»-^] , [•-•,(D-»] 



• • •], [ 



The first sets Ti r can be generated using a computer. The number of forests of the thirteen first sets Ti r is given in 
the following table: 



r 1 2 3 4 5 6 7 8 9 10 11 12 13 
card"H r 1 3 6 15 33 80 186 454 1109 2772 7012 18053 47066 



(119) 



We also define the set Ti of the forests made of trees belonging to Q of size larger or equal to 2, with edges relabeled 
by ii, «2i ■ • ■ so that all the labels are different. We call Ti r the subset of T-L made of the forests h such that \h\ — h = r. 
We note that for h £ H r , the number of edges in h is precisely equal to r, as in each tree of Q the number of edges is 
one less than the number of elementary nodes. On the contrary, for h £ T~L r , the number r is not equal to the number 
of edges of h since there is no edge inside the composite nodes. The two first sets T-L r are 




(120) 



where the . . . represent the forests obtained from the ones drawn by permuting the indices ij and by changing the 
directions of the edges. 

As in the case of the trees, we define for a forest h of T-L or T-L the set c(h) of the composite nodes of h. the set 
e(h) of the elementary nodes of h, and the set o(h) of the outer edges of h (i.e. the set of the edges of h linking 
two elementary nodes from different composite nodes, and not the edges linking two elementary nodes of the same 
composite node). 



B. Function W^'^ defined on the forests 

Most of the functions which were defined previously on the trees of Q and of Q in section IV Bl naturally extend 
to forests. We choose a map 8 from Ti to Ti which preserves the forest structure of the composite nodes. The map 
8 transforms a forest h G Ti r into a forest h! £ Ti r . It roots each tree of h, adds an internal tree structure to the 
elementary nodes inside each composite node, transforms the edges between composite nodes into edges between 
elementary nodes, and adds a direction and a label to all of the edges. 

The operator * applied on a forest h £ Ti is defined as the forest constituted of the trees of h to which the operator 
* has been applied to: the forest h* is thus made of trees having all their elementary nodes linked to the elementary 
root. 

The functions U v , v and V£ can now be extended naturally to the forests of Ti. For arbitrary numerical functions ip, 



24 



r\ and £, we define 



\e£e(h) J \eGe(?i) 

Vs(h)= n e(m(o)). 
The function also extends naturally to the forests. For h £ H r , we define 



(121) 
(122) 



w fw = E -E EW0W)^(W 



(123) 



For the same reason as in the case of the trees, we have another equivalent definition for applied to a forest: 



wf{h) = E - E twwm(W) 



iiGZ i r £2 



(124) 



As in the case of the trees, W^(h) does not depend on the choice of the map 9 if the function £ is even. 



C. Explicit formula for the cumulants of the current 



Equation (jlOlj) gives the generating function of the cumulants of the current £(7) in terms of a parameter B, which 
can be determined from the equation (|98p giving 7 in terms of B. This parametric representation of £(7) does not 
give a direct access to the cumulants of the current. To obtain them, the parameter B must be eliminated between 
equations (|101[) and (|98[) . This elimination can be done perturbatively in 7. At the first orders in 7, it gives explicit 
expressions for the first cumulants of the current. 

The use of the Lagrange inversion formula (see e.g. [42| or |43|) allows to go further, by inverting equation (|98|) 
explicitly to all orders in 7. This gives an expression of B as a formal series in 7. Inserting this expression of B in 
equation (|101[) finally gives £(7) as an explicit formal series in 7. The calculation is explained in appendix [Gj where 
we obtain the following result: 



E(l) - J7 1 - x ^ Y 



L - 1 ^ r! 

r=2 



r-1 



E 



Sf(h) 



heu r 



(125) 



The functions 77 and £\ were defined previously in equations and We observe that the sums over the trees in 
the parametric equations (|101[) and for £"(7) have became sums over forests. The forest symmetry factor Sf(h) 
is defined by 



r n 7> (N-i)i n |g| 



S/(/0 v p f (h) 11%) 



(126) 



The tree symmetry factor £((7) was defined in equation (|76|) . The factor Pf(h) is equal to the cardinal of the 
permutation group of the identical trees of h. If the forest h is made of mi times a first tree, mi times a second tree, 
. . . , then Pf(h) is simply equal to the product of the factorials of the mj. 
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Up to order 5 in 7, the expression (|125[) of £'(7) gives (using the notation W for W^f* 

(9W[*-*-m] - W[@] - 9W[ 



+ 



p(l - x) 
L 3 7 4 



48 



W 



72 



5W[ 



L 4 j 5 
28800 



900W[* 





+ 75W 


4 




[-<} 


• « 

4 


» • 

> 



450W [•-©-•] +27W[@] -3600W[< 

1200VK [(§)-•,•— •] - 2025W[«— •— •,« 
- 8100W[«-«-«, •—•,•—•] - 900W[(3), 



- 1200TU 

] +450VF [•—•—•,( 
3150IU 



-300VF[(3) • •] 



25IU 



+ G(7 6 ) ■ (127) 

We observe that this expression of E(j) allows to recover the known formulas (jlip and (fT2j) of the diffusion constant 
and of the third cumulant of the current. We also recover the expression (|19[) announced for the fourth cumulant 
E 4 {x). 



D. Thermodynamic limit with vanishing asymmetry 



For a system with particle density p = n/L fixed, it is possible to take the limit of the expression (|125[) of E(j) if 
the asymmetry 1 — x is of order 1 / y/~L. We set 



\0gXy/Lp(l - p) 



If 1 — x goes to zero when L goes to infinity, we have 

1-x - 



2 

2$ 



(128) 



(129) 



We want to replace x by its expression in terms of 3> in the function £\ used in (|125[) . From the definition of $, we 
can write 



A 



tanh 



if z = 
if z ^0 



We also want to take the limit L — > 00 in the function rj used in (|125|) . From Stirling's formula, we have 



r){z) 



1-P 



g 2p(l-p)i 



(130) 



(131) 



From the two previous equations, using the definition of WJ'^ff), the expression (|125[) of £^(7) becomes in the 



limit L 



r=2 



2 r-l r | 



E E- E 

h&V.T-i ti£Z tr-i£2 



[/ -2 (fl(M*) x V 

zi->z 2 ,2i->e 2 P (l-p)l, zi->tanh 



... \(*(>0) 



(132) 
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We set Uj = ij/y/Lp(l — p). When L goes to oo, the Riemann sums over the indices ij become integrals. We finally 
obtain 



E[l) - Jl ~{1-X)Y h r m(p(l - p))(r+m L 3{r-V)/2Y 



(133) 



with 



f iy-1 roo r oo U (9(h)*) X V z ^ tanh n z i^\(d(h)) 

M „ . tiL- ,„,... ^ — S/W ■ (-4) 

For r = 2 and r = 3, we recover the functions h r given in section Til CI after a few simplifications. The exponential 
e~ z / 2 in the definition of h r (G>) ensures the convergence at infinity of the integrals. The convergence at Uj = is 
more complicated to prove. We checked it explicitly for r between 2 and 5. For general r, it should be related to the 
cancellation of the arbitrary constant A from the finite size expressions. 

We note that the expression (|133[) of the cumulants of the current is in fact also true if 1/vT <C 1 — X, as all 
the expressions remain bounded and nonzero in the limit $ — > oo. It gives the cumulants of the current in the 
Kardar-Parisi-Zhang regime (see section Hi C[) as 

^ ~ (1 - x)h r (oo)(p(l - p ))(r+l)/2 L 3(r-l)/2 _ (135) 
P 

The case 1 — x <C l/V~L corresponding to the limit 4> — > is more complicated. From an exact calculation for r = 3 
[3ij], and from a numerical evaluation of the integrals in h r ($>) for r — 4, we conjecture that for r > 3, the limit $ — > 
of h r (&) gives 

KW - 2 T\ 2 ^ 1 » (136) 
(r-l)\ 

where the are the Bernoulli numbers. This gives for the r-th cumulant of the current (replacing $ by its expression 
in terms of 1 — x) 

E r - (^yyU - - p)) r L 2r - 2 . (137) 

We observe that this expression matches exactly the one obtained by taking the limit v — > oo in the expression of 
the cumulant E r in the weakly asymmetric scaling 1 — x ~ vjh [35|. We conjecture that the expression (|137|) gives 
the correct value of the r-th cumulant of the current in all the intermediate regime 1/L<C 1 — .t<C 1/vL, which lies 
between the Edwards- Wilkinson regime 1 — x <C 1/L and the Kardar-Parisi-Zhang regime — x. Proving 

that (|137p does hold in this whole regime would require the calculation of finite size corrections to the Riemann 
integrals. It was done in [34j in the case of the third cumulant. 



VIII. CONCLUSION 



The fluctuations of the stationary state current in the periodic asymmetric exclusion process have already been 
studied much in the past. The diffusion constant was first obtained in the totally asymmetric model [27| . and the 
result was then generalized to the case of the partially asymmetric model (28l . [33| . Higher cumulants of the current 
have also been calculated, both in the totally asymmetric model (30l |3"H and in the partially asymmetric model 




In this paper, we generalized these previous results for the cumulants of the current. Solving the functional Bcthe 
equation of the model with an arbitrary asymmetry between the hopping rates of the particles, we obtained exact 
finite size expressions for all the cumulants of the current. These expressions are given in terms of sums over tree sets, 
and have a nice combinatorial structure. It would be interesting to obtain them in a more direct way, starting from 
a graph theoretic formulation of the problem such as the one used in 0, EH to describe the stationary measure. In 
particular, having a combinatorial interpretation of the binomial coefficients, of the trees, and of the indices ij that 
appear in the expressions for the cumulants would be nice. Another interesting question is whether such combinatorial 
structures also appear in the case of the open ASEP and of the ASEP with several classes of particles. 
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Taking the thermodynamic limit in our exact formulas for the cumulants of the current in various scalings for the 
asymmetry, we observe three different regimes for the fluctuations of the current. In the regime of weakest asymmetry, 
which contains the model with symmetric hopping rates, it is known that the dynamics of the system is described at 
large scales by the Edwards- Wilkinson equation. Similarly, in the regime of strongest asymmetry which contains the 
totally asymmetric model, it is also known that the system is described at large scales by the Kardar-Parisi-Zhang 
equation. These two regimes are separated by an intermediate regime. It would be interesting to know how the 
system can be described at large scales in this regime. It would also be nice to know whether the presence of this 
intermediate regime between the Edwards- Wilkinson and the Kardar-Parisi-Zhang regimes is a general feature or is 
particular to the periodic exclusion process. 

The fluctuations of the current have also been studied for the exclusion process on the infinite line Z, using random 
matrix techniques. We emphasize that these results are not given by simply taking the infinite size limit in our 
expressions for the cumulants of the current, as the long time limit giving the stationary state does not commute with 
the large system size limit. Finding the crossover between the current fluctuations on the ring and on the infinite line 
is still an open question. 

We note that our calculation of the cumulants of the current involves finding the "ground state" of the non hermitian 
Hamiltonian of the XXZ spin chain with A > I and twisted boundary conditions. Our solution relies on a perturbative 
expansion of the twist parameter near a point where all the Bethe roots are equal. A generalization of this method 
to the calculation of other eigenvalues would be useful. 
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Appendix A: Proof of the functional equation "beyond the equator" 



In this appendix, we construct a polynomial P of degree L — n related by the functional equation Q37p to the 
polynomial Q of degree n solution of the functional Bethe equation (|32p . This construction is a straightforward 
generalization to the case of nonzero twist of the construction in the article [4l| by Pronko and Stroganov. 

For a system with n < L/2, we first divide the functional Bethe equation (|32p by Q(t/x)Q(t)Q(xt). We obtain 

R(t) _ e^(l-t) L , x n (l-xt) L 



Q(t/x)Q(xt) Q(t/x)Q(t) Q(t)Q(xt) ' 

If all the zeros of the polynomials Q(t/x) and Q(xt) in the variable t are different, which is generically the case, the 
rational function (I — t) L / (Q(t / x)Q(t)) can be written in the form 

^ - U ® ■ V ® , Wit)i (A2) 



Q(t/x)Q(t) Q(t/x) Q(t) 

where the polynomial W is of degree L — 2n, and the polynomials U and V are of degree at most n — 1 . In terms of 
the polynomials U, V and W, the functional equation (|AI [) reads 

R(t) _ e^Vit) + x^Uixt) , e^Uit) + ^) + + _ 



Q{t/x)Q{xt) Q(t) Q(t/x) Q{xt) 

If all the zeros of Q(t) arc different from the zeros of Q(t/x)Q(xt), which is again generically the case, the lhs of 
the previous equation has no pole in the zeros of Q{t), unlike the rhs. The previous equation then implies that the 
coefficients of these poles vanish. Since the polynomials U and V are of degree strictly lower than the degree of Q, 
we must necessarily have 

V(t) = ~x n e- L ~<U(xt) . (A4) 



With the help of the previous equation between the polynomials U and V, we can now eliminate V in equation (|A2|) . 
We obtain 

^ - U ® XnU ^ , w{t) _ (A5) 



Q(t/x)Q(t) Q(t/x) e L ~> Q{t) 
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It is now natural to write the polynomial W in the form 



W(t) = X(t) - ^X(xt) , (A6) 



with a polynomial X of degree L — 2n. If we call Wk the coefficient of t k in W(i), then the polynomial X is given by 

L -2™ ... 4-k 

*(*) - E r^k^ ■ (A7) 



k=0 

In terms of the polynomial X, equation (|A5|) becomes 

(l-t) L ( U(t) \ x 11 (U{xt) 



Q(t/x)Q(t) " y Q (t/x) ' Xit) ) ~ ^ XM) +X{xt) ) ' (A8) 
We now introduce the polynomial P, of degree L — n, defined by 

P(t) = (1 - x n e- L i)Q(0)(U(xt) + Q(t)X(xt)) . (A9) 
In terms of the polynomial P, equation (|A8j) finally becomes 

(1 - .T™ e - L7 )Q(0)(1 - t) L = Q(t)P(t/x) - x n e~ L ^Q(t/x)P(t) . (A10) 
We have recovered equation (|37j) . In particular, in t = 0, this equation implies 

P(0) = 1 . (All) 



Appendix B: Proof of the expression of w(t) as a sum over rooted trees 

In this appendix, we prove the expression (|68[) of w(t) as a sum over rooted trees. We start from the functional 
equation (|66|) and from the expansion in powers of C of w(t) (|67|) . The expansion of e rX ^ w ^ in powers of C gives 



e rXM*)] 



exp rY,X[w k {t)]C k = £ - £ XK(t)]C 

\ fc=l / s=0 ' ' \fc=l 



OO 



E^y E ■■• E ^KW].^kW]c ll+ - +t ' 

s=0 ' fei=l k 8 =l 

OO s , \ OO 

£7 e Ul... n^[».c)P)'' 



s=0 h.l 2 ,h....e 

h+h+h+:-s 

= f> •£ f|<K^. (B1) 

s=0 h,h,t 3 ,...€ti, i=l *" 
h+h+h + ---=s 

Between the second and the third line in the previous expression, we called U the number of indices kj such that 
kj = i. The multinomial coefficient ^ ) counts the number of s-tuples of indices (fci, . . . , k s ) corresponding to 
the same sequence of indices (Ji, 1%, h, ■ ■ ■)■ 
The equation (|66j) for w(t) now reads 

r=l s=0 ii,ia,ls,...eN, i=l *' 

(rodd) ; 1+ ; 2+ ; 3 _|_...= s 

At order k in C, this expression gives Wk(t) in terms of the Wi(t) for i < k. It means that the previous equation allows 
to solve recursively w(t) in powers of C. At first order in C, we find w(t) = h{t)C + O (C 2 ). We associate to the 
term at first order in C the tree © of Q° . More precisely, defining x(@) — h{t) and 5 (©) = 1, we have 



Mt) -sm~^m- (B3) 
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We will now show by recurrence on the integer k that it is possible to define suitable xid) an d S (g) for all the trees 
in so that 

-»(«)- E^- < B *> 

geg° ° VJJ 

The factor S (g) is a symmetry factor associated to the tree g, while the factor x(g) contain all the dependency on 
the asymmetry x of Wk{t). If we assume that for any i < k, Wi(t) is given by (|B4[) . then 

%(f 1 v / h \ ffi *[x(gf )] m ? 

2^ l m W m W 111 „ , r<i,™W 



m< i) + ...+mfj? 0| =Z i 



y TT *W )] . (B5) 

«' ...,m S„,6N J =1 m j ! > 3 



'IS. 



We called gj 1 ' 1 , . . . , 9\go^ the trees of the set Q°. The multinomial coefficient ( m w '* m (0 ) counts the number of terms 
in the expansion of X[wi(t)] li such that the trees g^ appear rrij times. We obtain from (|B2|) 



oo 



«*(*) = E E E E 

r+; 1 +2( 2 +3; 3 + ...=fe m f),..., m w eN 

I w 2 I 

mf'+...+mfJ ) ;| =i 1 
m< 2 »+...+mJ"> =i 2 



is? I , \ / oo ia°i 




v - n n „„ - : 1 1 w n n ^v? 1 i ■ m 



We will now associate a tree g G to each value of the tuple (r, s, Zi , 2 2 , '3, • ■ • , mf, • ■ • , TO [gj>p m \ 1 • ■ • 7 TO [g°| 1 ■ • ■) m 
the previous expression for Wk{t). The (odd) integer r will be the size of the composite root of g, that is, the number 
of elementary nodes contained in the same composite nodes as the elementary root (including the elementary root). 
The integer s will be the number of composite nodes directly linked by an edge to the composite root of g. The tree 
g is then made of a composite root of size r to which the trees <?) are attached rrij times for all i and j. With the 
notations of the previous equation, we observe that holds if we define x(d) an d S (g) by the recursive formulas 

x(g) = h(ty J] Xbc(g')}, (B7) 

g subtree of g 
connected to the 
composite root 



and 



1 _ (-l)^(H!)V TT TT m . 

Sob) r*r! llll m W !5o(5 « r « " ^ > 

As we observed on the expansion (|61l) of w(t) up to order 5 in C, x(ff) is obtained by labeling the composite nodes c of 
g by ft,(t)l c l, by labeling the edges of g by the X operator, and by reading the tree from the composite root, applying 
X on the term coming from the subtree attached to the corresponding edge, and multiplying the terms coming from 
different subtrees attached to the same composite node. 
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The recursive expression for S a (g) can be written explicitly as 



So(g) 




(-l)^(|c|!!) 2 | c | 



||\2| |l+number of subtrees attached to c in g 



|c|3|c|! 



n 



gi distinct subtrees 

attached to the 
composite root of g 
(with multiplicity mi) 



mi! 



n 



g2 distinct subtrees 

attached to the 
composite root of g\ 
(with multiplicity 1712) 



m 2 \ 



(B9) 



We see that we can rewrite it as 



So(g) 




-l)^(|c|!!) 2 |cr 



f size of the 
composite root 1 

^ of 3 



(number of equivalent^ 
choices of the ) 
composite root of g J 

Pis) 



(BIO) 



where P(g) is the number of permutations of the composite nodes of g leaving it invariant (except for the position of 
the root). The size of the composite root of g comes from the fact that the number of neighbors of the composite root 
is equal to the number of the subtrees attached to it in g, whereas the number of neighbors of any other composite 
node is equal to 1 plus the number of the subtrees attached to it in g. But, the number of equivalent choices of the 
composite root is equal to the number of equivalent choices of the elementary root of g divided by the size of the 
composite root of g. Thus, we can write 



1 

S4g) 




(-l)^(|c|!!) 2 |cr 



/number of equivalent 
choices of the 
elementary root of g 



{ 



P(g) 



(Bll) 



Appendix C: Generating function of the trees 

In section IIV1 it was shown that w r (t) is given by a sum over the trees of Q r with the symmetry factor S defined 
in equation (|76|). In this appendix, we study another interesting example of such a sum over Q r with the symmetry 
factor S: the "generating function of Q r " , which will be noted Z r (z). It is used in section IVIB I to recover the known 
expression for the current fluctuations in the totally asymmetric limit. It is defined by 

-card c(g) 

Z r ( Z ) = £ ~S(gT ' (C1) 
where card c(g) is the number of composite nodes of g. We will show that this generating function is given by 

Z r {z) = ^—T{[r{z + l)-2j] . (C2) 
r x r! J -- L 

In order to prove this expression, we first replace the sum over unrooted trees by a sum over rooted trees, by doing 
the reverse of the reasoning we used to go from the expression to the expression ([77)1 of Wk it) ■ We obtain 

1 z cardc(g') 

* w "'jSw- (C3) 

We now consider the generating function 

00 „card c(g ) /^card e{g ) 

U(C,z) =Y j rZ r {z)C r = E ^ . ( C4 ) 

r=i s'ee= bo[9 > 
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where carde(g') is the number of elementary nodes of g' . From a calculation similar to the one performed in appendix 
[B]to extract from equation (|59")l the expression (|68p of w(t), the generating function U = U(C, z) is the solution of 



U = z arcsinh(Ce' 7 ) 



The parameter C can be written in terms of U 



C = e- u sum 



U 



The Lagrange inversion formula (see e.g. [42|,|43|) then gives U as 



u(c, Z )=y:- 



r=l 



sinh r (u/;z) 



(u-i) 



cr= \yL_ 

zitt ^-^ r 

r=l 



(111 



2r e ru(z+l) 1 2 



Co (e 2M/z - If ' 



(C5) 



(C6) 



(C7) 



where the contour Co surrounds but not the other poles of \j(e 2u l z — 1). We perform the change of variables 
w = e 2u l z in the contour integral. The contour of integration C\ now surrounds 1. Wc find 



smh r (u/z) 



r(z + l) 

zdw 2 r w 2 j 



2i7r id 2^ - l) r 



2 

r - 1 



- 1 



Wc finally find 



00 r—l 



U(C,z) = z^2 — l[(r(z + l)-2j) 



(C8) 



(C9) 



r=l j = l 



We recover the expression (|C2|) for Z r (z). We checked this expression for r between 1 and 16 by constructing 149388 
trees of Q and by computing their symmetry factors. 

Appendix D: Proof of the equivalence between the two definitions of WJp' 

In this appendix, we show that both definitions (|89|) and (|90|) of W^(g) are equivalent. This will be illustrated on 
several examples of trees g £ Q . 

If g is a tree from Q r , W%'^(g) can be expressed as a sum over all integer values (both negative and nonnegative) of 

the indices i\, . . . , i r -i of the product of the functions U^^ and V£ applied to the trees 9(g) and 0(g)* elements of Q r . 
Taking 8(g) as the argument of U v>v and 9(g)* as the argument of corresponds to the definition (|59"|) of W^(g), 
while the opposite choice corresponds to the definition ([90)) of W%'*(cf). Starting from (|90|) . changes of variables can be 
performed on the indices i\, . . . , i r -i- In general, these changes of variables are such that U v . v (9(g)*) is transformed 
into an expression which can not be expressed as U v ^ v (g') for some tree g' £ Q r . The same is true for V^(9(g)). 
However, there exists changes of variables corresponding to transforming the tree 9(g)* into a tree gi £ Q r and the 
tree 9(g) into a tree g[ £ Q r as we will show in the following. We will see that it is in fact possible to find a sequence 
of such changes of variables of the indices ij corresponding to a sequence of doublets of trees (gj , g'^ ) beginning with 
(9(g)*, 9(g)) and ending with (9(g), 9(g)*): 



Wf(g) = E ■ • ' E U v , n (9(9r)V ( (9(g)) (definition flM} of Wf(g)) 

iiGZ i, — i£Z 

= E ■ • ■ E u <pMvM = ■ ■ • = E ■ • ■ E u <p,Mv*(sD 



E ■ • ■ E U <P*Q(9))VdW) (definition © of Wf(g)) 



(Dl) 
(D2) 
(D3) 



This will prove the equivalence between the two definitions l{89]) and (|9Q]l of 
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1. Changes of variables corresponding to changes of trees 

We consider the expression 

]T... ]T EW<7i)n(si) , (D4) 

i, — i£Z 

where g\ and g[ are trees from Q r . The changes of variables ij — > — ij simply correspond to transforming the trees g± 
and g[ into g 2 and g' 2 obtained from g\ and g[ by reversing the direction of the edge labeled by ij. On the contrary, 
other changes of variables of the ij, only correspond to changing U VyTI (gi) into U v _ n (g 2 ) in the previous expression, but 
do not necessarily correspond to transforming V^g^) into an expression which can be put in the form V^(g 2 ). Other 
changes of variables can correspond to changing V^(g' 1 ) into V^{g' 2 ) but do not transform U v ^(gi) into an expression 
of the form U VyTI {g 2 ). In this subsection, we will study separately the changes of variables changing g\ to g 2 (while 
transforming V^(g[) in an arbitrary way), and the ones changing g[ to g' 2 (while transforming U v ^(gi) arbitrarily). 

We begin with the changes of variables modifying the tree g\. If an elementary node e of g\ is surrounded by 
k edges pointing to e (labeled by a\, . . . , a k ) and by I edges pointing away from e (labeled by 61, . . . , b{), then e 
contributes to U v ^ n (gi) only through £(e), which is equal to the sum of the aj minus the sum of the bj\ 

\ai 61/ 

"' • )\ ' '" => e(e) = {a 1 + ... + a k )-(b 1 + ... + b l ) . (D5) 

a y 

Changing one of the into minus the sum of some of the aj (different from ai) plus the sum of some of the bj is 
equivalent to removing these aj and these bj from £(e) and to add them to £{e'), where e' is the elementary node which 
is linked by the edge labeled to the elementary node e. Thus, it is equivalent to detach from e the edges labeled by 
the relevant aj and bj and to attach them back to e'. This change of variables corresponds then to a transfer in the 
tree g\ of some edges from e to e'. This transfer is local in the sense that the nodes e and e' are neighbors. Similarly, 
the change of variables transforming one of the bi into 6j minus the sum of some of the bj (different from bi) plus the 
sum of some of the aj corresponds to a transfer of edges from e to the node e" such that the edge between e and e" 
is labeled bi in g\ . 

We now continue with the changes of variables modifying the tree g[. We call a\, . . . , a k the indices labeling the 
edges which point away from the elementary root of g[ and b\, . . . , bi the other indices, which point to the elementary 
root. Let p be an edge (inner or outer) of g[ labeled by one of the a^. The contribution of a; to V^(g) comes from all 
the m(o) for o an outer edge of g[ located between the composite root of g[ and the edge p (including m{p) if p is an 
outer edge): 

v d9i) = £( fl 4 + a 5 + a 6 - bi - b 2 ) 

x£(a 2 + a 4 + a 5 + a 6 - b x - b 2 ) ■ (D6) 
x (terms that do not contain a$) 

We consider the change of variables changing to at minus the the sum of the aj (different from a^) labeling the 
edges of the subtree of g[ beginning at the edge p plus the sum of the bj labeling the edges of the same subtree. This 
change of variables consists in detaching from g[ all the subtrees beginning at the neighboring edges of p farther from 
the elementary root than p, and attaching them to the elementary root of g[ . Contrary to the changes of variables 
considered previously for g\ , this change of variables is non-local: we attach to the elementary root of g[ some subtrees 
initially arbitrarily far from it. Similarly, we can consider an edge p of g[ labeled by one of the bi. The corresponding 
change of variables consists in changing bi into bi minus the sum of the bj (different from bi) labeling the edges of the 
subtree of g[ beginning with the edge p plus the sum of the aj labeling the edges of the same subtree. The meaning 
of this change of variables for the tree g[ is the same as for the previous change of variables involving bi instead of a^. 

Depending on the respective structure of the trees g\ and g[ , it is possible that the two sets of changes of variables 
that we introduced in the two previous paragraphs have a non-empty intersection. These changes of variables then 
preserve the structure of the expression (|D4|) . modifying both g\ and g[. These are the changes of variables that we 
will use in the following to show that both definitions of are equivalent. 
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2. Example 1: linear tree with composite nodes of size 1 

, and by making for 6(g) the 



We begin by considering the case where g is the linear tree g 
following choice: 

Q 



6(9) 



t5 
> 

•u 



is and 6(g)* 




(D7) 



The definition ((90j) of W^(g) gives 



Wf(g) = J2---J2 U ^n 




x Vz 



(D8) 



Vi'V 

■ • t P T l{h,i2,i3,ii,is,—i5 — U — h -h — i\) 

iiSZ i 5 GZ 

x €(ii,h + *2j *i + ii + «3, *i + h + iz + Uj «i + «2 + «3 + U + is) 



where we used the notations 



ipr](ai, ...,ai) = (y>(oi) + • • • + <f(a>i)) x 77(01) x . . . x 77(0;) 
£(&i,...,& m ) =^(61) x ... x £(b m ) , 



(D9) 
(D10) 



in order to lighten the expressions. In the previous expression of W%'^(g), we will make the following sequence of 
changes of variables: 



7-5 -ll - 12 - 73 - *4 + 75 



7 4 — S- -7i - «2 - 7 3 + 7 4 , 



(Dll) 



73 — » — i\ - 7 2 + 73 , and i 2 -> — ii + 7 2 
The first change of variables 75 — > —i\ —72 — 73 — 74 + 75 gives 

= ^2 ■ --^2 ( P r ?(*1^2^3,i4,7 5 - 7 4 - 7 3 - 7 2 - 7'l,-7 5 ) 



X + 7 2 , 7i + 7 2 + 7 3 , 7'i + 7 2 + 7 3 + 7 4 , 7 5 ) 



ilGZ i 5 GZ 




1> 



'2 



(D12) 
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The second change of variables i± — > —i\ — ii — 13 + 14 leads to 

W^(g) = ^ ■ ■ .y] frj(i 1 ,i 2 ,i3,i4 - is ~ 12 - h>h ~ U, — *5)^(*i»*i + *2, «i + «2 + «3,*4,«5) 




(D13) 



By the third change of variables 13 — > — ii — 12 + 13, we obtain 

W^(g) = ^2 . . 'PV(ii,h,i3 ~ h ~ h, U - hJ5 - U, -is) x + ii-, *3i H, is) 

iiGZ i 5 ez 




(D14) 



Finally, the fourth change of variables 12 —ii + «2 gives 

^^(5) = y] • ■ ^ W(*l,*2 - *1,*3 - «2,«4 - *3) *5 - «4, -£5) X i 2 , l 3 , «4, £5) 




(D15) 



which is the expression corresponding to the definition (|89[) of W%'*(jg). 

We observe that at each step, we recover the "upper part" of the tree g (starting at i§) in the tree on the left, and 
the "lower part" of the tree g (starting from i\) in the tree on the right. At each step, we destroy the uppermost part 
of the tree on the right and we rebuild it in the tree on the left: the complete structure of g is thus preserved between 
the two trees. 

The previous proof for the tree g can be easily extended to a linear tree of arbitrary size. The introduction of 
composite nodes of size strictly larger than 1 does not change much if 9(g) is chosen such that the tree structure of 
the elementary nodes is linear. 



3. Example 2: branched tree with composite nodes of size > 1 

We now consider the example of the tree 




(D16) 
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and we choose 9(g) such that 




and 6(g)* 




(D17) 



Starting from 




(D18) 



[*5/ V6 

Ha 

= <pT}(ii,i2,ia, ii, — «S, *6,«7, -«8, — «i — «2 — «3 - *4 + «5 - *6 - *7 + is) 

£(*3, «7, -«8, «4 ~ «5 + «6 + *7 ~ *s) , 

we will step by step destroy the upper part of the tree on the right while rebuilding it in the tree on the left, similarly 
to what we did in the previous section for the case of the linear tree. We begin with the change of variables i\ — > ii — 13, 
12 — > ii — ii + is — is — ii + is- We obtain 

W^{g) = (pr)(ii - i 3 , i 2 - it + h - is - «7 + is, h, U, -is, h, h, -is, ~h - 12) 

x £(«3, ii, ~is, *4 - is + H + ij - is) 




E-E^ 



The change of variables 24 — > 14 + 15 — is — ii + is then leads to 

W^(g) = <pn(ii - i 3 , «2 - U, 13, ii + ib - is - ii + is, -is, is, 17, -is, ~h - «2)C(«3, «7, -is, ii) 

( \ / /T> \ 



(D19) 




E-E^ 



Finally, we perform the change of variables i§ — > 15 + 17, is — > ie + is and we obtain 

w ^ (9) = <PV(h,ii ~ h, -h ~ h, «2 - U, ii + is - is, -is - %7, is + h, h, -h)^(h,U, %7, -is) 



(D20) 



E-E^. 





(D21) 



We observe that for any tree, the procedure that we have presented to transform the expression ([50]) into the expression 
([89]) of (<?) works, which proves the equivalence between these two definitions. We note that we do not need at 
this stage to have a function £ even. This is only necessary to ensure that W^(g) does not depend of the map 9 that 
has been chosen (c./. appendix [E]) . 
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Appendix E: Independence of W^'^ with respect to the choice of the map 9 

We prove in this appendix that the function defined in (|89|) and (|90f is independent of the choice of the map 

6 provided that the function £ is even. We recall that the map 6 transforms the trees of Q into trees of Q. We must 
show that for any tree g £ G, W^(g) is independent of the choice of the labeling and of the directions of the edges 
of 6(g), of the choice of the elementary root in 6(g), and of the choice of the internal tree structure of the composite 
nodes of 6(g). We will see that each of these properties can be easily seen from either one of the two definition (|89|) 
and (|90p of W^(g), the equivalence of which is proved in appendix [Pi 

We begin with the independence of W^(g) with respect of the choice of the labeling of 6(g). Permuting the 
indices which label the edges of 6(g) in ([89| or (|90|) is equivalent to permuting the indices ij in U ip ,, 1 (6(g))V^(6(g)*) 
or Utp tT] (d(g)*)V^ (6(g)). The summation over the ij implies then that W^(g) stays unchanged by the permutation 
of the labels of the edges. 

We continue with the independence with respect to the choice of the directions of the edges. Reversing the directions 
of the edge labeled by ij is equivalent to replacing ij by — ij in U V}TI (d(g))V^(6(g)*) or U v _ v (6(g)*)V^(6(g)). Again, 
the summation over the ij in the definition of W^(g) ensures the independence with respect to the choice of the 
directions of the edges in 6(g). 

We now move on to the independence with respect to the position of the elementary root of 6(g). It is now useful 
to consider the definition of WR'^(g) instead of the other one. By definition, U v ^(6(g)) does not depend on the 
position of the elementary root in 6(g). On the contrary, V^(6(g)*) depends on it: depending on the position of the 
root of 6(g), the orientation of some of the edges of 6(g)* can be modified. This corresponds to changing some of the 
ij which appear in V^(6(g)*) into — ij. But this transformation of some ij into — ij only takes place in V^(6(g)*) and 
not in U ip ^ v (6(g)). Thus, we can not use the summation over the ij to conclude. We need to use the following remark: 
the factor V^(6(g)*) depends on the indices ij only through £,(ij) (there is no linear combination of the indices ij when 
is applied on the tree 6(g)*). The transformations of V^(6(g)*) resulting from the modification of the position of 
the root of 6(g) are thus equivalent to modifying some factors £,(ij) into £(— ij). It is thus crucial that the function £ 
is even so that the independence of WR'^(g) with respect to the position of the root in 6(g) is verified. 

We still must show the independence of WR *(g) with respect to the choice of the internal tree structure of the 
composite nodes by the map 6. We will use here the definition (|90|) of W^(g). By definition, V^(6(g)) does not 
depend on the internal tree structure on the composite nodes of g (but depends however on the orientation of the 
edges in this internal tree structure). The tree 6(g)* does not depend either on the internal tree structure of the 
composite nodes of g. We thus conclude that WR'^ does not depend on the choice of the composite nodes of 6(g). 

We have finally proved that W^(g) does not depend on the choice of the map 6. The fact that the function £ is 
even is only needed to prove the independence with respect to the choice of the root of 6(g). 

Appendix F: Proof of the parametric expression of the polynomial Q 

In this appendix, we perform the expansion of the expression (|77[) of w(t) in powers of t, which leads to the 
expression of the polynomial Q(t) in terms of the function defined in (|89p and ([517)1 . 

1. Explicit expression for w(t) 

We start with the expression ([77)1 of w(t): 

H *w = EiEx(ft), (fi) 

g£Gk Vy; j=l 

where the function \ is defined in equation We recall that the trees gj <G G% are the trees obtained from g by 

choosing for elementary root the j-th elementary node of g (for some ordering of the elementary nodes of g). We will 
write w(t) in a more explicit form by expanding all the functions h(t) = (1 — t) L /t n in powers of t, except the one 
corresponding to the elementary root. 

We first label the k elementary nodes of g € Qk by ii, . . . , ik so that in the tree gj, the index ij labels the elementary 
root. Then, for all m different from j, we expand the h(t) corresponding to the m-th elementary node as 



h(t) = hij 



(F2) 
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For g € Gk, we have 

k k 

Exfe) = M*)£ E ^■••Vi^ +I ---M <1+ "-" H, - 1 " Hj+1 "'" H *^(ffi)» (F3) 

j=l j=l ti,...,»i_i,*i + i * fc €Z 

where V(gj), which corresponds to the action of the operators X in x{9j)i IS equal to a product of factors of the form 
£\(z) with z a sum of indices from i\, . . . , %k and 

/- / i A if z = 

^) = (i±g if ^o ' (F4) 

More precisely, V(gj) is equal to the product over all the composite nodes c of gj (except the composite root) of £a 
applied to the sum of the indices labeling the subtree of gj made of the composite nodes that can be attained from c 
by moving away from the composite root. For example, for the labeled tree 

14 / J 3 * 




9 = 



(F5) 




«6- 



we have 

V(gi) = V(g 2 ) = V(g 3 ) = 6(i 4 )k(i 5 + i e + i 7 + is + isO&ftO&fo) (F6) 

V{9i) = £a(«i + «2 + i 3 + «5 + k + h + h + «9)£a(*5 + k + h + is + i9)£,\{is)£\(i9) (F7) 

V(gs) = V(g 6 ) = V(g-) = &(i x + z 2 + z 3 + ^Ka^Ka^KaM (F8) 

V{gs) = £a(«i + «2 + «3 + ? : 4 + h + H + ii + £9)6^9)6^1 + «2 + «3 + ? ; 4)£a(«4) (F9) 

^(59) = + «2 + «3 + «4 + «5 + *6 + *7 + »8)£\(i8)£\(il + «2 + h + h)£,\{h) ■ (F10) 

We now rename by ij in equation (|F3|) . We obtain 

k k 

Exfe) = M*)E E ^-^.li^-^-'WjiUi) ■ (in' 

j=l j=l ii,...,i fc _iGZ 

We will now use the following property (which can be seen easily on the example of the labeled tree g defined in (|F5I) , 
using the fact that the function £ A is even): 

Vigjh^i, = Vigk)^-^-...-^ . (F12) 
Then, for j ^ k, we perform the change of variables ij — > —i\ — ... — ik-i- We find 

k fc-1 

E^fe) =/l (*)E E ^1 ■ ■ ■h ij _ 1 h-i 1 -...-i h _ 1 h ij+1 .■■h ik _ 1 t-**V(g k ) 

3 = 1 3 = 1 ii,...,i fc _iGZ 

E (F13) 

„ ii ••• ife — 1 

!l,...,!t_l€* 



Factoring the product of the h m , we can write 

Exfe)=M<) E 



1 



h - pi /1 - fifc-i /; +— ii — — ifc_i 



ft il ...ft jfc _ 1 A- il -...- <fc _ 1 V(fffc) • (F14) 



Finally, we find for u;(i) 



"'A' 



5(o) ^ \h i A il "' ft,-, ,i ifc -! h-i,- *-n---t*-i 



x ^ ...fc ik _ 1 /i_ il _..._i k _ 1 V(fffc) • (F15) 

We emphasize that the choice of the particular rooted version of g that we used here (i.e. gk) is completely arbitrary: 
any rooting of g could have been chosen. 
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2. Explicit expression for a(t) 



In order to calculate E(pf), we need to extract the negative powers in t of w(t) to find a(t). Thus, we only need to 
keep the negative powers of h(t)/t z = (1 — t) L /t n+z . We will use the notation [■ • • ]( \ for the negative powers in t of 
a sum of powers of t. We have 



h(t) 



n+2— 1 

= E 

-) i=0 



(F16) 



It is convenient for the calculation of the generating function of the cumulants of the current to write this as a scries 
in (1 — t) instead of a series in t. We have 



hit) 



(-) 



n-\-z — l oo 

= E E 

i=o ;=o 



L 



% — n — z 
I 



(-i) i+l (i-ty 



(F17) 



The sum over i can be performed in the previous expression by considering the term of order in u in the product of 
the two formal series 



/=() 

OO 



1 =y 

(1 



This gives the binomial sum 



n-\-z— 1 

E 

i=0 



L\ fi — n — z 



-3 - 1 
/ 



(-l)«u 



and we find 



Using this and 



h{t) 



= E 

(-) ;=o 



(_1)*+J 



L-Z- 1 

n + z — 1 



L-Z- 1 
n + z — 1 



(-1) 



n-\-z— 1 



n + 2 



(-1) 



n-\-z 



we have finally found the following expression for a(t): 



oo oo 1 / ( L—l — X ' 

fe=0 ;=0 geGk W ii,...,i h _ieZ \ U+iJ 



/£-l-l\ 

f L ) 

Vn+i fc _i/ 



L 

n + ii 



L 



n + ik—ij \n — H 



( L ^ ) 

\n—l—i 1 —... — i k - 1 J 

( . L . ) 
\n—%i—...—%k-i/ 

L 



(F18) 
(F19) 



(F20) 



(F21) 



(F22) 



V(ff fc ) 



(F23) 



The global factor 2 comes from the factor 2 in the definition ([55]) of w(t). 

We note that the previous equation can be expressed in terms of the tree functions defined in section IV Bl Indeed, 
V(gk) is equal to V^ x (9(g)) if the map 9 from Q to Q of section IV Bl is defined such that the choice of the elementary 
root in #(17) is the same as in g^, and if all the arrows are pointing away from the elementary root. With such a map 
9, the product of binomial coefficients in the previous expression is also equal to U i p ltV (9(g)*), the functions tpi and 77 
being defined in equations (|9Tj) and (|92|) . Thus, we can write for a(t) the following expression 



00 00 

ait) = 2j2J2i-i) kn c k ii-t) 1 J2sT) E u Vl ,Mgy)VtMg)) 

geg k ii,...,i fe _iez 



(F24) 



fc=0 1=0 
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Introducing the function defined in (|90[) and using the expression (|62j) of a(t) in terms of the polynomial Q, we 

finally obtain 



where B is given in terms of the parameter C defined in (|58|) by 



S = 2{-l) n ( L \c . (F26) 



Appendix G: Proof of the explicit expression of the generating function E(y) 



In this appendix, we prove the explicit expression (|125|) for the generating function of the cumulants of the current 
E(j). We show that this expression, which involves a sum over forest sets, is a consequence of the parametric 
expression (|98|) and (|10ip of E(j) involving sums over tree sets. 



1. Elimination of B in the parametric expression of E(-y) 



The equations (|98|) and (|101[) express respectively 7 and E(j) as formal series in a parameter B. These equations 
are of the form 



f(B) = J2f k B\ 



k=l 



and 



Eh) - Jl 



P 



= g(B) = J29kB k 



(Gl) 



(G2) 



k=2 



We want to invert the relation (|G1[) to obtain an expression for B in terms of 7, and then insert this expression 
into (|G2|) to obtain an explicit expression in 7 of E(j). The inversion of the relation (jGlj) between i? and 7 can be 
performed using the Lagrange inversion formula (see e.g. f42l. |43|). This formula implies that, for a formal series g 
such as the one defined in (|G2[) , we have 



[<?(£)] (7-) = - 



/(5) r 



(G3) 



(B-i) 



where [• • • ]( 7 rj and [• • • stand respectively for the coefficient of the term in 7 r and the one of the term in B 1 

of the expression • • • . Using the expansion (|G2[) of g(B) into powers of B, we obtain 



From ([G2]), we find for £(7) 
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(G4) 
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We need to expand f(B) r into powers of B for r strictly positive integer. By the multinomial theorem (with a 
negative exponent), we have 
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The term in 5~ fc of j{B)~ r is then 
1 
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(B- k ) 
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and we find for £'(7) (shifting k of 1 
E{l)-Jl = f-I (2 



r=2 



1 ( fjl 



Ei i E^ 1 ^ E (r-i + fc+fc + .-oiIIn - 



fc=i 
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3=1 
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We must now replace in the previous equation the coefficients fj and gu+i by their respective values given by ([98 
and poip . We have in particular f\ = —1/L. We finally find, putting together all the powers of 2: 
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We have obtained an expression of E(j) in terms of tuples of trees. We now want to replace these tuples of trees by 
forests. 



2. Expression of £(7) as a sum over forests 

Equation (|G9|) involves a sum over b\ trees of size 2, . . . , bj~-i trees of size fc, 1 + trees of size k + 1, 6/c+i trees 
of size fc + 2, ... Wc note / the m-tuple (ordered set of m trees) constituted by these trees, ordered by increasing 
size. We have m = Y^jLiftj + ^j,k)- In equation (|G9[) . we thus sum over all the m-tuples of trees / such that the 61 
first trees of / are of size 2, the 62 following of size 3, . . . The sum over the m-tuples / can be replaced by a sum over 
the forests h containing b\ trees of size 2, . . . , bk-i trees of size k, 1 + bk trees of size k + 1, bk+i trees of size k + 2, 
. . . There, we must take into account the fact that several m-tuples of trees will correspond to the same forest: two 
m-tuples correspond to the same forest if and only if they are made of the same trees but not in the same order. We 
consider the action on the set of the m-tuples of trees / of the subgroup of the group of permutations of m elements 
keeping the size of the trees in increasing order in the m-tuplcs. The number of m-tuples of trees corresponding to 
the same forest as a given m-tuple / is equal to the cardinal of the orbit of / under this group action. Thus, we have 



II ('', + hk) 



card Orb(/) - 3 1 - , (G10) 

where the product of the factorials of the bj + Sj^k is the cardinal of the group of the permutations leaving the sizes 
of the trees in increasing order in the m-tuple /, while Pf(h), defined after equation (|126p , is equal to the cardinal of 
the stabilizer of / under the group action. Calling go one of the trees of size k + 1, we obtain 

(6 fc + l)(r- I + 61 + 63 + ...)' 
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We now want to put W^^ 2 (go) together with the product for g of size k + 1 (g ^ go) of the W^\(g). We use the 
fact that, by the definition (|124|) of acting on a forest, we have 

v w^mwZ^m W^\{g bh ) 
h r ^ S(g ) S( 9l ) S(g bk ) 

h={g ,...,g b AeU 
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|ffo|=---=|sb fe I 



(^■■■^-.(n-W*)', E -ssr--^ 



(k + l) b »W^ z2 ({g ,...,g bk }) 
b k + l S(g )x ...xS(g bk ) ' 



(G12) 



with the notation {go, . . . , g bk ] for the forest made of the trees go, ■ ■ ■ , g bk ■ The factor {k + l) bk is the product of the 
factors k + 1 of all the W^'^ (ffj), while the factor l/(bfc + 1) comes from the symmetrization of z i— > z 2 over all the 
trees of the forest {50 , ■ • ■ , 9b k }- We obtain 

E{l)-Jl _ 2{l-x) ^ 1 / L 7V'^ (fc + l) bfc+1 x (r -! + &!+ b 2 + ...)! 

^ V ; r=2 v 7 fc=l 61,62,— GN I ( — I) 3 
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geh 

* k m us.r 1 — ■ ™ 

Vj bj + 5i fe trees " 
of size j + 1 

We now change 6fc + 1 into so that for all j, we have exactly bj trees of size j + 1 (even for j = fc). Then, we can 
put the sum over k inside the sum over the forests. We can write 

E(j) - Jj 2(1 - x) ^ 1 / L 7 \ r x - (r-2 + h + b 2 + ...)! 
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The number of trees in the forest {g <G h, \g\ = k + 1} being equal to b k , W TI ^ 2 ({<? G ^7 Iff I = k + 1}) is equal to zero if 
bk = 0. The constraint l6 fc >o can thus be forgotten. By definition of for the trees and the forests, we have then 



of size j + 1 

J> + l) b * x W££ a ({$ e ft, \g\ = k+ 1}) W^i(g) = ( + ^ ] ^^({.9 € fe» • (G15) 
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(G16) 



For j k, the factor (j + 1)^ comes from the W^\{g) for g one of the bj trees of size j + 1. We have then 

E{l) -Jj _ 
P 

2(i-») ^i / DyV v v 

of size j + 1 

We observe that r — 1 = &i + 2&2 + ... is equal to \h\ — h and that b\ + &2 + • • • is equal to the number ft of trees in 
ft. It implies that (r — 2 + &i + 62 + . . .)! is equal to (|ft| — 1)!. We also note that the j + 1 are the sizes of the trees 
g E ft. We can write 

We have recovered the forest symmetry factor of h defined in (|126[) . We can remove the sum over the bj and sum 
directly over the forests h € H r -i- We finally obtain the explicit expression (|125[) for E(-f): 

m-J 1 = 2(1 - X ) {-iyuY v w^ 2 (h) 

p L(L-l)2-t 2 r r\ ^ S f (h) ' 1 ' 
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